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(57) A method for allowing a number of wireless transceivers to exchange information (data, voice or video) with each 
other. A first frame of information is multiplexed over a number of wideband frequency bands at a first transceiver, 
and the information transmitted to a second transceiver. The information is received and processed at the second 
transceiver. The information is differentially encoded using phase shift keying. In addition, after a pre-selected time 
interval, the first transceiver may transmit again. During the preselected time interval, the second transceiver may 
exchange information with another transceiver in a time duplex fashion. The processing of the signal at the second 
transceiver may include estimating the phase differential of the transmitted signal and pre-distorting the transmitted 
signal. A transceiver includes an encoder for encoding information, a wideband frequency division multiplexer for 
multiplexing the information onto wideband frequency voice channels, and a local oscillator for upconverting the 
multiplexed information. The apparatus may include a processor for applying a Fourier transform to the multiplexed 
information to bring the information into the time domain for transmission. 
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replacing conventional analog devices. The system 
operates with relatively small cells. In other 
aspects, dynamic channel allocation and voice 
activation may be used to improve the capacity of the 
5 system. 

Advantages of the present invention include: 

1. It can be used indoors as well as outdoors using 
the same transceivers. If data is to be exchanged, as 
opposed to voice, the transceiver preferably contains 

10 an estimator to allow pre-distortion and post- 
distortion of the transmitted signal. 

2. The system, as compared with prior art systems 
omits the clock or carrier recovery, automatic gain 
control or passband limiter, power amplifier, an 

15 equalizer or an interleaver-deinter leaver, and 
therefore has low complexity. 

3. The system offers good speech quality, as well as 
low probabilities of dropped and blocked calls. It is 
robust against Doppler and multipath shifts. It is 

20 also robust against both impulse noise and narrowband 
interference. 

4. The system is flexible, such that at the expense 
of increased complexity of the DSP receiver it can be 
applied over noncontiguous bands. This is accomplished 

25 by dividing a 100 MHz (in one of the exemplary 
embodiments described here) band into several subbands 
each accommodating an integer number of voice 
channels. 

5. The system offers low frame delay (less than 26.2 
30 ms in the exemplary cellular embodiment described 

here). The transceiver requires low average 
transmitted power ( of the order of 20pW in the 
exemplary cellular embodiment described here) which 
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means power saving as well as enhanced biological 
safety. 

6. The system offers up to a 38 fold increase in 
capacity over the North American Advanced Mobile Phone 
5 System (AMPS) which uses analog frequency modulation. 

Operation of the system in accordance with 
the techniques described in this disclosure may permit 
compliance with technical requirements for spread 
spectrum systems. 

10 There is therefore disclosed in one aspect 

of the invention a method for allowing a number of 
wireless transceivers to exchange information (data, 
voice or video) with each other. In the method, a 
first frame of information is multiplexed over a 

15 number of frequency bands at a first transceiver, and 
the information transmitted to a second transceiver. 
In a cellular implementation, the second transceiver 
may be a base station with capacity to exchange 
information with several other transceivers. The 

20 information is received and processed at the second 
transceiver. The frequency bands are selected to 
occupy a wideband and are preferably contiguous, with 
the information being differentially encoded using 
phase shift keying. 

25 A signal may then be sent from the second 

transceiver to the first transceiver and de-processed 
at the first transceiver. In addition, after a pre- 
selected time interval, the first transceiver 
transmits again. During the preselected time interval, 

30 the second transceiver may exchange information with 
another transceiver in a time duplex fashion. 

The processing of the signal at the second 
transceiver may include estimating the phase 
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differential of the transmitted signal and pre- 
distorting the transmitted signal. 

The time intervals used by the transceivers 
may be assigned so that a plurality of time intervals 
5 are made available to the first transceiver for each 

time interval made available to the second transceiver 
while the first transceiver is transmitting, and for 
a plurality of time intervals to be made available to 
the second transceiver for each time interval made 

10 available to the first transceiver otherwise. 
Frequencies may also be borrowed by one base station 
from an adjacent base station. Thus if one base 
station has available a first set of frequencies, and 
another a second set of distinct frequencies, then a 

15 portion of the frequencies in the first set may be 
temporarily re-assigned to the second base station. 

In an implementation of the invention for a 
local area network, each transceiver may be made 
identical except for its address. 

20 Apparatus for carrying out the method of the 

invention is also described here. The basic apparatus 
is a transceiver which will include an encoder for 
encoding information, a wideband frequency division 
multiplexer for multiplexing the information onto 

25 wideband frequency voice channels, and a local 
oscillator for upconverting the multiplexed 
information. The apparatus may include a processor for 
applying a Fourier transform to the multiplexed 
information to bring the information into the time 

30 domain for transmission. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

There will now be described a preferred 
embodiment of the invention, with reference to the 
drawings, by way of illustration, in which like 
numerals denote like elements and in which: 

Figures la and lb are schematics of a prior 
art receiver and transmitter respectively; 

Figure 2 is a schematic showing the use of 
the available frequencies according to one aspect of 
the invention for use with cellular applications; 

Figure 3a is a schematic showing an 
idealized pulse for transmission over a cellular 
system; 

Figure 3b is a schematic showing a modified 
version of the pulse shown in Figure 3a; 

Figure 3c is a schematic showing a further 
modified version of the pulse shown in Figure 3a; 

Figure 4 is a schematic showing an exemplary 
protocol for cellular communication; 

Figure 5a is a block diagram showing the 
structure and function of an embodiment of the 
transmitter of a cellular portable in accordance with 
the invention; 

Figure 5b is a block diagram showing the 
structure and function of an embodiment of the 
transmitter and receiver of a cellular base station in 
accordance with the invention; 

Figure 5c is a block diagram showing the 
structure and function of an embodiment of the 
receiver of a cellular portable in accordance with the 
invention; 

Figure 6a is a flow diagram showing the 
function of the processor in either of Figures 5a or 
5b; 
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Figure 6b is a schematic showing the 
function of the deprocessor in either of Figures 5b or 
5c; 

Figure 6c is a schematic further 
illustrating the operation of the processor and 
deprocessor shown in Figures 6a and 6b; 

Figure 7a is a schematic showing the 
structure and function of the channel estimator in 
Figure 5b; 

Figure 7b is a flow chart showing the 
operation of the channel estimator of Figures 5b and 
7a; 

Figures 8a , 8b and 8c are respectively 
schematics of 126, 63 and 7 cell reuse patterns; 

Figures 9a and 9b are schematics showing 
transmit protocols according to one aspect of the 
invention; 

Figure 10 is a schematic showing the use of 
the available frequencies according to another aspect 
of the invention for use with local area network 
applications; 

Figure 11a is a schematic showing an 
idealized pulse for transmission over a local network 
system; 

Figure lib is a schematic showing a modified 
version of the pulse shown in Figure 11a; 

Figure 11c is a schematic showing a further 
modified version of the pulse shown in Figure 11a; 

Figure 12 is a schematic showing a preferred 
protocol for local area network communication; 

Figure 13a is a block diagram showing the 
structure and function of an embodiment of the 
transmitter of a local area network transceiver 
according to the invention; 
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Figure 13b is a block diagram showing the 
structure and function of an embodiment of a further 
local area network transceiver according to the 
invention; 

5 Figure 13c is a block diagram showing the 

structure and function of an embodiment of the 
receiver of a local area network transceiver according 
to the invention; 

Figure 14a is a flow diagram showing the 
10 function of the processor in either of Figures 13a or 
13b; 

Figure 14b is a schematic showing the 
function of the deprocessor in either of Figures 13b 
or 13c; 

15 Figure 14c is a schematic further 

illustrating the operation of the processor and 
deprocessor shown in Figures 14a and 14b; and 

Figure 15 is a schematic showing the 
structure and function of the channel estimator in 

20 Figure 13b. 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

Introduction 

The benefits of the invention can be readily 
appreciated with reference to Figure 1, which shows a 

25 prior art transmitter/receiver configuration for a 
portable unit. The transmitter includes a vocoder 110, 
an interleaver 112, a modulator 114, a filter 116, 
local oscillator 118, power amplifier (PA) 120 and 
antenna 122. The receiver includes an LNA 124, a local 

30 oscillator 126, a filter 128, automatic gain control 
(AGC) 130 with an associated passband hardlimiter not 
separately shown, carrier recovery 132, sampler 134, 
clock recovery 136, adaptive (or fixed) equalizer 138, 
demodulator 140, deinterleaver 142 and decoder 146. 
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With implementation of the present invention, several 
of the blocks shown in Figure 1 are not required. 
These are the interleaver 112, deinterleaver 142, 
power amplifier 120, automatic gain control 130 with 
passband hard-limiter , clock recovery 136 and carrier 
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recovery 132, and the equalizer 138. It will now be 
explained how the proposed system obtains the omission 
of these blocks without impairing the quality and 
capacity of the system, 
5 in this disclosure there will be described 

two systems as examples of the implementation of the 
invention. The system described first here will apply 
to a cellular system with a number of portable 
transceivers and base stations (BS). Then will be 

10 described a local area network implementation. A local 
area network will typically be a system of equal 
transceivers. The invention may also be implemented 
with combinations of cellular and local area networks, 
or to a system with a number of equal transceivers and 

15 a master or controlling transceiver. "Equal" as used 
here means that the transceivers have more or less the 
same processing equipment and processing capabilities. 
The system described here is primarily for the k 
exchange of voice information. 

20 Link set-up and termination protocols 

between transceivers, and the equipment required to 
implement them, are well understood in the art as well 
as the basic structure of radio transceivers that may 
be used to implement the invention. Hence these 

25 elements are not described here. What is described 
here are the novel operational, functional and 
structural elements that constitute the invention. 

Cellular Implementation of Wideband Modulation 

30 The present invention proposes in one 

embodiment a wideband modulation scheme for exchange 
of information between transceivers such as portables 
and base stations* 
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Wideband in this patent document is 
described in the context of Wideband-Orthogonal 
Frequency Domain Modulation (W-OFDM or wideband OFDM) - 
In OFDM, the entire available bandwidth B is divided 
5 into a number of points K, where adjacent points are 
separated by a frequency band Af , that is B - KAf . The 
K points are grouped into a frame of K x points and two 
tail slots of K 2 points each, so that K = 1^ + 2K 2 . The 
frame carries the information intended for 

10 transmission under the form of multilevel differential 
phase shift keying (MDPSK) symbols or differential 
quadrature amplitude modulated (DQAM) symbols. Thus 
each point in the frame corresponds to one information 
symbol. The two tail slots act as guard bands to 

15 ensure that the out-of-band signal is below a certain 
power level. For example, when a pul6e P(f) is 
selected for pulse shaping and the out-of-band signal 
has to be ydB or less relative to the in-band signal, 
K 2 is selected such that 

20 20.1og 10 |P(f )/P(0) | s; y for f * K 2 Af . 

When the pulse is a raised-cosine pulse with a roll- 
off p and when the number of levels each symbol can 
take is M, the bit rate is equal to I^logjM/^t + (1 
+ P)/Af) where (1 + p)/Af is the duration of the frame 

25 and 6t is the guard time required to take into account 
the delay of arrival and the delay spread due to 
raultipath. In this case, the bandwidth efficiency, 
which is defined as the ratio between the bit rate and 
the bandwidth, is equal to: 

30 log 2 M/((l + P + 6tAf)(l + 2^/^)) 

In wideband-OFDM, both K and Af are selected 
sufficiently large to achieve a high throughput as 
well as to reduce the effects on the BER of the clock 
error, the Doppler shift and the frequency offset 
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between the LO in the transmitter and the one in the 
receiver. To show what is meant by "K and Af are 
selected sufficiently large" , consider the effect of 
increasing K and Af on (1) the clock error, (2) the 
5 Doppler shift and (3) the frequency offset between the 
LO in the transmitter and the LO in the receiver. 

(1) When a clock error at a transceiver of 
value z occurs in the time domain, it causes a shift 
in the phase difference between adjacent symbols in 

10 the frequency domain of value 2nAfr. When t is equal 
to %T where T is duration of one time domain sample 
and % is any real value, the shift is equal to 2rrAfxT. 
Hence, r causes a shift in the phase difference 
between adjacent symbols of value 2u%/K 1 since T is 

15 equal to 1/{K X L±). By doubling the number of symbols 
from K± to 2K a the shift in the phase difference is 
reduced by half from 271%/^ to 71%/^. Thus, the effect 
of the clock error on the BER is reduced by increasing 
K. 

20 (2) When there is relative motion between 

the transmitting transceiver and the receiving 
transceiver, a Doppler shift occurs with a maximum 
absolute value |V/A| where V is the relative velocity 
between the two transceivers and X is the wavelength 

25 of the travelling wave corresponding to the carrier 
frequency f c (i.e. f c is the frequency corresponding 
to the middle point in the frame). Such a Doppler 
shift causes a sampling error in the frequency domain 
of the same amount, or equivalent ly, it causes a 

30 sampling error of V/(AAf) relative to one symbol 
sample. Thus, the effect of the Doppler shift on the 
BER is reduced by increasing Af . 

(3) When a frequency offset between the LO 
in the transmitter and the one in the receiver occurs 
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with a value f 0 , it causes a sampling error in the 
frequency domain of the same amount, or eguivalently, 
it causes a sampling error of f 0 /Af relative to one 
symbol sample. Thus, the effect on the BER of the 
5 frequency offset between the LO in transmitter and the 
one in the receiver is reduced by increasing Af . 

In summary, OFDM with a K and a Af large 
enough to be able to achieve a specific throughput and 
large enough to be able to avoid using either a clock 

10 or a carrier recovery device without substantially 
affecting the BER is referred to here as Wideband- 
OFDM. As an example, let us assume that MDPSK is used 
in an OFDM system with the number M of levels, with a 
carrier frequency f c , with a raised cosine pulse of 

15 roll-off p, with the LO at the receiver having a 
frequency offset f Q relative to the LO at the 
transmitter (so that the frequency offset between the 
carrier frequencies in the first and second 
transceivers of the multiplexed information is f Q ), 

20 with a given maximum expected clock error r = %T at 
the receiving transceiver, where T is the duration of 
one time domain sample, and with a maximum expected 
relative velocity V between the transceivers. Thus, 
in order to ensure that the out-of-band signal is ydB 

25 or less relative to the in-band signal and to be able 
to avoid using either a clock or a carrier recovery 
device without substantially affecting the BER we have 
to: 

1. Find the acceptable sampling error Af ' , relative 
30 to one symbol sample, which does not 

substantially affect the BER. This can be done 
using the following rules: 

When 0.2 £ P < 0.3, Af = 7.50% 
When 0.3 * 0 £ 0.4, Af = 10.0% 
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When 0.4 s 0 s 0.5, Af = 12.5% 

When 0.5 <; p s 0.6, Af = 15.0% 
2. Find Af such that: 

V/(AAf) + f Q /Af £ Af 
5 3. Find K 2 such that 

20.1og 10 |P<f )/P(0) | s y for f * K 2 of 
4. Find K x such that 

2n%/K 1 < tf/M 
In this case, we refer to OFDM as Wideband-OFDM. 
10 Element 4 is a necessary condition for wideband OFDM, 
and given a sampling error, the sampling error may be 
corrected with the methods described in this patent 
document . 

To implement wideband modulation, Orthogonal 

15 Frequency Division Multiplexing (OFDM) is preferred in 
which the information, for example encoded speech, is 
multiplexed over a number of contiguous frequency 
bands. Wideband OFDM forces the channel to be 
frequency selective and causes two types of linear 

20 distortion: amplitude distortion and phase distortion. 

To reduce the effect of amplitude distortion the 
modulation is preferably phase modulation, while the 
effect of phase distortion is reduced by employing 
differential phase modulation. Hence the modulation 

25 may be referred to as Differential OFDM (DOFDM). 
Unlike in other proposed schemes, neither pilot tones 
nor diversity are required in DOFDM . Possibly, 
quadrature amplitude modulation might be used, but 
amplitude modulation makes it difficult to equalize 

30 the distorting effects of the channel on the signal. 

To implement wideband modulation in a 
cellular system with a plurality of portables and one 
or more base stations, a 100MHz band is divided into 
4096 points, as shown in Figure 2, plus two tail slots 



2064975 



15 

of 195.3 KHz each. The 4096 points represent N voice 
channels (vc). Adjacent points are separated by 24.414 
KHz and each point represents a Differential eight 
Phase Shift Keying (D8PSK) Symbol e^< n >, where C(n) = 
5 C(n-l) + 0(n) + X( n )« ^( n ) takes one of the eight 
values {0, 2n78, 4jt/8, 14tt/8} with equal 

probability for n = 1, 2, . .., 4096 and 0(0) takes an 
arbitrary value. %(n) also take6 an arbitrary value. 
%(n) may be used as a security key and will be known 

10 only to the transmitter and receiver. Information in 
the form of output bits from a vocoder are mapped onto 
#(n) . Vocoders are well known in the art and do not 
need to be described in detail here. The focus here is 
to transmit the bits with an acceptable Bit Error 

15 Rate, i.e. with a BER £10" 2 for voice and <10~ 8 for 
data. 

Ideally, 40,96 p& (=1/24.414 KHz ) is the 
minimum duration required for one frame to be 
transmitted without frequency domain intersymbol 

20 interference. This can be achieved using a Raised 
Cosine (RC) pulse with zero roll-off, as shown in 
Figure 3a. Figure 3a illustrates a rectangular (time 
domain) window corresponding to the RC (frequency 
domain) pulse. Such a pulse, however, requires an 

25 infinite frequency band. To alleviate such a 
requirement, an RC pulse with a 20% roll-off (i.e. 0 
- 0.2) may be used as shown in Figure 3b. The frame 
duration has increased by 20% to 49.152 ps. The two 
tail slots of 195.3 KHz each (i.e. 8 points each) 

30 ensure that the signal outside the entire band of 
100.39 MHz is below -50 dB. To allow the frame to 
spread over the time as a consequence of the multipath 
nature of the channel, an excess frame duration of 
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2.848 ps is provided as shown in Figure 3c, making the 
frame duration 52 ye in total. 

Since the frame duration is 52 ys, the frame 
rate is 252 frames per 13.104 ms or equivalently, 126 
5 full duplex frames may be transmitted/received every 

13.104 ms. The reason for pre-selecting an interval of 
13.104 ms is to ensure a transmission delay to allow 
one transceiver to communicate with other transceivers 
at the same time, but must not be so long that the 

10 delay becomes unacceptable to the user. Delays longer 
than about 40 ms are too great for voice, and it is 
preferable to be lower. For data, the delay may be 
longer and still be acceptable. 

In the exemplary embodiment described here, 

15 three bit rates are considered for the vocoder: 18.77 
Kbps, 9.16 Kbps and 6.18 Kbps. Table I displays the 
structure of a vc slot and the number N of vc for each 
vocoder rate. The control symbols in each vc slot are 
required for handoff and power control. Figure 2 shows 

20 that N vc can be transmitted simultaneously. This is 
known as Frequency Division Multiple Access. Figure 3c 
shows that 126 full duplex frames can be transmitted 
every 13.104 ms in a Time Division Multiple Access 
fashion (TDMA). The total number of Full Duplex voice 

25 channels (FDvc) is therefore 126xN and is shown in 
Table I. 

To ensure that the channel is slowly fading, 
a Time Division Duplex protocol for exchange of 
information between the portable and the base station 
30 is proposed as illustrated in Figure 4. The protocol 
is as follows: 

1. The portable transmits a frame 410 over one vc 
slot. See the discussion in relation to Figure 5a 
below. 
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2. The Base Station (BS) receives the frame 410 from 
the portable and processes (analyzes) it as shown and 
discussed in relation to Figure 5b below. 

3. Based on the received signal, the BS predistorts 
5 a frame 420 and transmits it to the portable over the 

same vc slot, 520 ps or some other suitable time 
interval later in which the channel does not change 
substantially. The time interval will depend on 
factors such as the frequency,, speed of the 
10 transceiver and other environmental factors. 

4. The portable receives the frame from the BS. See 
the discussion in relation to Figure 5c below. 

5. Steps 1 through 4 are repeated, as for example by 
the transmission of the next frame 430, every 13.104 

15 ms until the call is terminated. 

During 520 ps, a portable travelling outdoor 
at 100 km/hr moves 1.44 cm, which leaves the outdoor 
channel largely unchanged. Indoors, a portable moving 
at 2 m/s moves 0.1 cm again leaving the channel 

20 unchanged. Assuming that the channel is reciprocal and 
stationary over 520 ps, a predistorted signal, 
transmitted by the BS, should reach the portable 
undistorted. 

From Figure 4, one can see that the portable 

25 transmits/receives one FDvc every 13.104 ms, while the 
BS can transmit /receive up to 21 frames or 
equivalently up to 21xN FDvc every 13.104 ms. The 
frames 440 labelled frame 2... frame 21 are frames that 
may be transmitted to other portables. This implies 

30 that while one BS is processing its data over 520 ps, 
six other BS can communicate to their corresponding 
portables in a Time Division Multiple Access (TDMA) 
fashion using the same frequency bands. Also, during 
the 13.104 ms, or such other preselected time interval 
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that is suitable, the BS may communicate with one or 
more other portables. 

When a portable is stationary during a call, 
it is possible with high probability to have the 
5 transmitted signal centered with several deep 
(frequency domain) nulls, hence, causing speech 
degradation. Also, narrowband interference over the vc 
slot can deteriorate the speech. In order to avoid 
both situations, the signal is preferably frequency 

10 hopped into a new vc slot within the same (frequency 
domain) frame. This frequency hopping is ordered by 
the BS which is constantly monitoring the channel 
frequency response. Monitoring techniques, as well as 
frequency hopping, are known in the art, and not 

15 described here further. When an unacceptable speech 
degradation is first noticed by the BS a probation 
period is initiated and maintained for at least 10 
cycles (i.e. 10x13.104 ms) unless speech degradation 
has ceased. In other words, the probation period is 

20 terminated if speech degradation has ceased. Frequency 
hopping is then ordered at the end of the probation 
period. The period of 10 cycles is long enough to 
indicate the portable stationarity and is short enough 
to allow speech interpolation between unacceptable 

25 speech frames, hence maintaining good speech quality. 

As known in the art, the BS ensures that no collisions 
take place between hopping portables. 

Digital Signal Processing 

30 The transmitter/receiver block diagrams 

corresponding to the protocol in Figure 4 are shown in 
Figures 5a, 5b and 5c. Figure 5a corresponds to step 
1 in the protocol described above. Speech is provided 
to a vocoder 510 where the speech is digitized and 
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coded to create bits of information. The bits are 
provided to the modulator 512 which turns them into 
D8PSK symbols, with three bits per symbol. The D8PSK 
symbols are then processed in the processor 514 which 
5 is described in more detail in Figure 6a. The output 

from the processor is then filtered in low pass filter 
516/ upconverted to RF frequencies using local 
oscillator 518 and transmitted by antenna 520. Figure 
5b corresponds to steps 2 and 3. 

10 In Figure 5b, the received signal at the 

base station is filtered in a bandpass filter 522, and 
down converted by mixing with the output of a local 
oscillator 524. The average power of the downcoverted 
signal is monitored by a power controller 525 that 

15 adjusts the average power to the specifications 
required by the sampler 526. The adjusted 
downconverted bits are then sampled in sampler 526 to 
produce bits of information. The bits are then 
processed in the deprocessor 528, described in more 

20 detail in Figure 6b. An estimate of the phase 
differential of the received signal is taken in the 
channel estimator 530, as described in more detail in 
relation to Figure 7a and 7b below, and the estimated 
phase differential is supplied to a decoder- 

25 demodulator 532 to correct the received bits. The 
estimated phase differential is also supplied to a 
pre-distorter 534 in the transmitter. At the 
transmitter in the Base Station, the same blocks are 
incorporated as in the portable transmitter except 

30 that a pre-distorter is used to alter the phase of the 
D8PSK symbols to make the channel appear Gaussian 
(ideal) as opposed to a fading channel. The pre- 
distorter 534 receives a signal corresponding to the 
estimated phase differential of the channel. On the 
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(believed reasonable) assumption that the channel is 
reciprocal, the signal being transmitted is pre- 
distorted with the estimated phase differential so 
that the received signal at the portable with which 
5 the BS is communicating will be corrected for any 
phase distortion over the channel. The advantage of 
rendering the channel Gaussian is a large saving in 
the power required to achieve an acceptable BER. The 
initial power control 525 also sends a signal to the 

10 pre-distorter 534 to adjust the transmitted power to 
an appropriate signal level for the sampler 526 in the 
portable's receiver depending on the average power of 
the received signal. Thus if the average power is too 
low, the transmitted power is increased and if the 

15 average power is too high, the transmitted power is 
decreased. The power controller 25 may also be used in 
frequency hopping to monitor the average power of the 
received frequency and determine when frequency 
hopping need take place. 

20 Figure 5c corresponds to step 4, and shows 

the receiver of the portable, which is the same as the 
receiver in the BS except it does not include an 
estimator or a power controller. These are not 
required in the portable on the assumption that the BS 

25 will carry out the phase estimation and the power 
control. However, if desired, the portable may include 
these functions. 

Figures 6a, 6b and 6c illustrate the 
function and structure of the processor and the 

30 deprocessor respectively in the transmitter and 
receiver. Software for modelling the function of the 
processor in a general purpose computer is attached 
hereto as Appendix A, and for modelling the function 
of the deprocessor as Appendix B. 
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Figure 6a shows that the processor is a DSP 
implementation of an RC pulse shaping filter with a 
20% roll-off, followed by an inverse Fourier 
transform. The processor first inverse Fourier 
transforms the 4096 D8PSK modulated symbols output 
from the modulator. The transformed symbols are then 
triplicated as a group so that the total number of 
samples is tripled, with three consecutive groups each 
consisting of the 4096 transformed symbols. The 
triplication of the signal is illustrated in Figure 
6c, where the symbols are shown as first delayed and 
added together . Next, as shown in Figures 6a and 6c, 
the three groups are windowed by a Raised Cosine 
window with a roll-off of 0.2 centered in the middle 
of the three groups. In other words, the processor 
takes D8PSK symbols in, pulse shapes them and inverse 
Fourier transforms them. On the other hand, the 
deprocessor undoes what the processor did, i.e. it 
removes the pulse shaping, then Fourier transforms the 
received signal to obtain the original D8PSK symbols. 
The first two blocks in Figure 6b are similar to the 
second two blocks in Figure 6a except for two 
differences. The two differences are as follows. In 
the first block of the deprocessor, the repeated 
groups of symbols are partially overlapped as shown in 
Figure 6c on the right hand side. In the second block, 
a rectangular window is used instead of the Raised 
Cosine. In the preferred implementation, the blocks 
are repeated three times but other numbers of 
repetition may be used. 

Figures 6a, 6b and 6c show that the DSP 
blocks used in the processor are identical to the ones 
used in the deprocessor, except for a small change in 
the two transforms and a small change in the shapes of 
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the two windows. Thus the same hardware can be used by 
both the processor and the deprocessor. 

Figure 7a shows a block diagram of an 
example of a preferred channel estimator, and Figure 
5 7b is a flow chart showing the operation of the phase 
estimator. Each of the steps is carried out in a 
computing means that may be a special purpose computer 
or a general purpose computer programmed to carry out 
the digital signal processing described here, as for 

10 example with the software that is attached hereto as 
appendix C. Other methods of estimating the channel 
may be used that obtain an estimate of the channel 
group delay or phase differential of the transmitted 
symbols. However, a preferred implementation is. 

15 described here. 

The first block in Figure 7a estimates the 
envelope A(n) for n = 1, 4096 of the (frequency 

domain) samples transmitted over the fading channel as 
output from the deprocessor. The estimate A 1 (n) is 

20 the square-root of t^he sum of the squares of the 
quadrature (Q) and inphase (I) samples output from the 
deprocessor. The second block performs the operation: 
Aln(A»(t)) = (A'(n) - A' (n-l))/A» (n), for n = 2, 
4096, where A'(n) is the estimate of A(n). The third 

25 block performs a Hilbert transform operation 
H[Aln(A» (t) ) ] on the result of the second block. 
H[ Aln(A' (t ) ) ] is an estimate of lAw(n)l for n = 2, 
4096, where Ao(n) is the phase differential of 
the transmitted signal (« is the phase of the signal). 

30 The Hilbert transform is preferably carried out by 
taking the discrete fast Fourier transform of the data 
record, multiplying the positive frequency spectrum of 
the transform by -i (square root -1), and the negative 
frequency spectrum of the transform by i, and taking 



2064975 



23 

the inverse discrete fast Fourier transform. The 
result is a set of symbols representing an estimate of 
the phase differential of the received signal, as 
determined from it6 sampled amplitude envelope. 
5 Instead of a Hilbert transform, a different 

estimation may be made to estimate the phase 
differential. In this case, firstly, after the 
electromagnetic signal has been sampled, a series of 
data frames of a number of consecutive amplitude 

10 samples (A(t)) of the electromagnetic signal are 
constructed. These data frames are then segmented into 
segments [t lf t 2 ], where the amplitude of the 
electromagnetic signal is at least a predetermined 
number of dB less than its running mean, for example, 

15 lOdB. The following calculation is then applied to 

these segments of the amplitude samples: 

Att(t) -l/t 0 -1 ,. 

1 + <f /t 0 >* 

20 where t 1 = t - t min , t mln is the time in [t lf t 2 ] when 
A(t) reaches its minimum, t is the time from the 
beginning of the Begment, and t 0 is the time from the 
instant the amplitude of the electromagnetic signal 
reaches its minimum during the segment until the 

25 amplitude reaches double its minimum during the 
segment. In other words, the phase differential may be 
calculated from 

Aci(t) « -t 0 /(t 0 2 + t» 2 ). 
The polarity of Aw(n) is extracted using the 

30 last block shown in Figure 7a, The estimate so 
calculated does not provide the sign of the 
differential. This may be determined by known 
techniques, for example by adding the phase 
differential to and subtracting the phase differential 
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from the received phase (tan" 1 (Q/I)) and taking the 
sign to be positive if the addition results in the 
smaller Euclidean distance to the expected value and 
negative if the subtraction results in the smaller 
Euclidean distance to the expected value. 

Equivalently, for each sample n, the ideal 
phase closest to w( n )+A(o(n ) is determined and labelled 
o + (n), and the ideal phase closest to o)(n)-Ad>(n) is 
determined and labelled 6>_(n). The two sums P = 
s|o + (n) - {o(n)+A(i)(n) } | and N = 2|a> + - {&>(n)-A(i>(n) } | 
are calculated. If P < N, then u>(n)+A»(n) is used to 
correct the signal, and if not then ci>(n)-Ati>(n) is used 
to correct the signal . 

For simplicity of the estimator, the 
determination of the sign need only be carried out for 
phase differentials greater than a predetermined 
threshold. This will be in the vicinity of a fade and 
may be accomplished by segmenting the data record into 
a segment in which the phase differential is larger 
than a selected threshold and setting the remainder of 
the data record to zero. This computation may be 
carried out with a simple discrimination circuit or 
equivalent computing means in the estimator. 

The bias 6u> of the channel group delay is 
estimated by averaging Ata'(n) over n for n = 1, . 
4096 where Ao)'(n) is the measured value of Aa>(n). The 
estimates A 1 (n) and Aw'fn) are used directly in the 
predistortion filter in Figure 5b, while the estimates 
Au(n) and 6g> of the unbiased channel group delay and 
of the bias of the channel group delay respectively 
are used in the demodulator. 

The complexity of the processor-deprocessor- 
channel estimator is displayed in Table II. Complexity 
is measured in Mega Instructions Per Second (MIPS) 
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where one instruction is defined as one complex 
addition , one complex multiplication and a storage of 
one complex number. It does not include overhead. 

The complexity of the processor-deprocessor- 
5 channel estimator in the BS is computed from the 
complexity of the Inverse Fast Fourier Transform 
( IFFT) /Fast Fourier Transform (FFT)/Hilbert Transform. 
The complexity is 4096 x 12 x 4 x 21/13. 104ms for the 
BS. For the portable, it is computed from the 

10 complexity of the FFT/IFFT per vc: (32 x 5 + 64 + 128 

+ 256 + 512 + 1024 + 2048 + 4096 )2/13 . 104ras for the 
portable with a 6.18 Kbps vocoder. Such a complexity 
assumes that the A/D converter operates at 100 MHz 
with 12 bit precision. As seen in . Table II, the 

15 portable has smaller complexity due to the fact that 
the portable transmits/receives one vc in 13.104 ms 
and the BS transmits/receives up to 21xN vc in 13.104 
ms. 

20 Reducing Analog Complexity 

Comparing Figure 1 (prior art) and Figure 5, 
it will be seen that several conventional blocks are 
not used in the present invention, namely the 
inter leaver-deinter leaver, the Power Amplifier (PA), 

25 both the clock and the carrier recovery, both the AGC 
with its associated Passband hard limiter, as well as 
the equalizer. 

From the BS point of view, the interleaves 
deinterleaver is not required since the signal is 

30 predistorted before transmission forcing the received 
samples to be independent. From the portable point of 
view, the interleaver-deinterleaver is not required as 
a separate entity from the vocoder due to the fact 
that the channel is highly frequency selective, hence 
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the interleaving/deinterleaving can be applied 
implicitly in the vocoder over one vc, without a need 
for a separate time domain interleaver/deinterleaver . 
This eliminates excess speech delays associated with 
5 interleaving/deinterleaving between frames. 

The PA is not required since the cells can 
have, as shown later, a radius of up to at least 250 
m outdoors and 30 m indoors, if the transmitted power 
is up to 6 dBra. Such a power can be generated by the 

10 Local Oscillator (LO) without a need for a PA. It is 
important to avoid using a PA since DOFDM generates a 
time domain signal with non constant envelope. A power 
efficient class C PA cannot be used without distorting 
the signal. A class A PA can be used at the expense of 

15 power efficiency. 

A clock recovery device is not required 
since a sampling error in the time domain is 
equivalent to a phase shift in the frequency domain. 
The phase shift is a linear function of frequency. It 

20 contributes to the bias in the channel group delay. 

Such a bias can be easily estimated and removed as 
mentioned previously by averaging w'fn) over n in the 
frequency domain. Such an estimate is accurate as long 
as the sampling error is less than 0.2 ps or 

25 equivalently less than 20 samples (since in this case, 
the corresponding phase shift is less than nr), and as 
long as the number of points in one vc is large enough 
as it is here. 

A carrier recovery device is not required 

30 since a carrier offset in the time domain is 
equivalent to a sampling error in the frequency 
domain. For the chosen RC pulse, a sampling error of 
up to 10% of the duration of one pulse is acceptable. 
This implies that a frequency offset of up to 2.414 
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KHz is acceptable regardless whether it is due to 
carrier offset as low as 1 part in a million, i.e. as 
low as 1 KHz per 1 GHz. When a carrier frequency 
higher than 2.414 GHz is required, one can decrease in 
5 Figure 2 the number of points per 100 MHz or one can 
use an RC pulse with a rolloff larger than 20%. 

Neither an AGC nor a Passband hard- limit er 
are required since the level of the received power may 
be controlled constantly. This is achieved as follows: 

10 The portable transmits a frame. The BS receives the 
frame and predistorts a frame intended for 
transmission accordingly, assuming that the channel is 
reciprocal and stationary over 520 ps. This includes 
controlling the transmitted power according to the 

15 received power. The BS transmits the predistorted 
frame and simultaneously orders the portable to 
control its power. The order is conveyed using the 
control symbol in the vc slot (See table I). The 
degree of power control may be determined using the 

20 power controller 525; and the instruction for the 
inclusion of a power control symbol in the vc may be 
sent from the power controller 525 to the pre- 
distorter 534. 

One advantage of wideband modulation over 

25 narrowband modulation is that the wideband signal does 
not experience short term fading the same way the 
narrowband one does. The wideband signal is mainly 
affected by shadowing and other long term effects 
which vary slowly and are easily monitored from one 

30 frame to the other as long as the same vc slot is used 
by the portable to transmit and receive (i.e. as long 
as TDD is employed). 

Finally, conventional equalization, whether 
it is linear or nonlinear, is not required simply 
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because there is little or no ISI. Also, from the 
portable point of view, each received vc is 
predistorted by the BS. Hence, the channel can be 
modeled approximately as an ideal memoryless Additive 
5 White Noise Gaussian (AWGN) channel, assuming channel 
reciprocity and stationarity over 520 ys. From the BS 
point of view, since the received signal is not 
predistorted by the portable prior transmission, the 
channel estimator is used to reduce the effect of the 
10 channel group delay. 

Smaller cells 

As mentioned previously, the LO generates a 
6 dBm average power, hence the signal power 

15 transmitted by the BS over one vc slot is (6 dBm - 
101og 10 N dB) while the signal power transmitted by the 
portable over one vc slot is 0 dBm. Also, since the 
noise power over a 100 MHz band is -94 dBm, it is (-94 
dBm - 101og 10 N dB) over one vc. A typical noise figure 

20 at the receiver is 7 dB. The penalty for not using a 
matched filter in the receiver is 1 dB. Combining 
together the above figures provides the portable with 
an (86 dB - path loss in dB) received signal to noise 
ratio (SNR), while it provides the BS with an (86 dB 

25 + 101og 10 N dB - path loss in dB) received SNR. 

For a path loss of 75 dB, the radius of the 
urban cell can be 250 m while it can be 30 m for the 
indoor cell. Such a path loss provides the portable 
with a 17 dB received SNR, while it provides the BS 

30 with an (11 dB + 101og 10 N dB) received SNR. From the 
portable point of view, the channel can be modeled 
approximately as an ideal AWGN channel, hence the 17 
dB received SNR results in a 2x1 0" 3 BER. On the other 
hand, the channel can be pessimistically modeled as a 
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Rayleigh fading channel from the BS point of view. The 
corresponding BER are displayed in Table III which 
shows that the achieved BER is 1 4x1 0~ 3 . A BER < 10" 2 
is acceptable for voice. 

5 

Cell Pattern Reuse 

From Table I, the number of Full Duplex 
voice channels (FDvc) that can be transmitted/received 
per frame is 136 over 100 MHz , for a 6.18 Kbps 

10 vocoder. If the bandwidth is halved to 50 MHz, the 
number of FDvc per frame is reduced to 68, the noise 
floor is reduced by 3 dB and the number of full duplex 
frames that a BS can transmit/receive is doubled to 
42, leaving the frame duration, the number of frames 

15 per 13.104 ms and the processor /deprocessor complexity 
unchanged. 

Reducing the available bandwidth directly 
affects the cell pattern reuse. This can be explained 
as follows, assuming that we are required to offer a 

20 minimum of 136 FDvc per cell, that the vocoder rate is 
6.18 Kbps and that the cell radius is fixed at 1 km 
outdoors and 30 m indoors. For a 100 MHz band, we 
assign one frame per cell and offer 136 FDvc per cell. 
In this case, the cell pattern reuse consists of 126 

25 cells as shown in Figure 8a which displays a seven 
. layer structure. For a 50 MHz band, we assign two 
frames per cell and offer 136 FDvc per cell, hence 
reducing our cell pattern reuse to a 63 cell pattern 
as shown in Figure 8b which displays a five layer 

30 structure. If the available bandwidth is as low a6 
5.86 MHz, we have 8 vc per frame. Hence we have to 
assign 18 frames per cell in order to offer the 
minimum required number of FDvc per cell. This reduces 
the cell pattern reuse to as low as a 7 cell pattern 
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as shown in Figure 8c which displays a two layer 
structure. 

In Figures 8a, b and c, a shaded area is 
shown around the center of the pattern, indicating 19, 
5 38 and 126 full duplex frames that the central BS can 
transmit /receive respectively. Tables IVa, b and c 
show the number of cell layers in each cell pattern 
reuse, the coverage area in Km 2 of the pattern reuse 
for both the indoor and the urban environments, as 
10 well as the carrier to interference ratio (CIR) in dB, 
for the 100 MHz, 50 MHz and 5.96 MHz bands, 
respectively* In all cases, the CIR is large enough to 
sustain a toll quality speech. 

15 Transmission/Reception Protocol 

Since the number of FDvc a portable can 
transmit/receive is one, while the number of FDvc a BS 
can transmit/receive is much larger as shown in Table 
V for each of the three vocoder rates, we have chosen 
20 the following transmission/reception protocol: 

1. The portable transmits a frame over a vc. 

2 . Seven adjacent BS receive the frame from the 
portable. 

3. One BS transmits to the portable, depending for 
25 example on the strength of the received signal by each 

of the BS. 

The control of this protocol may use any of 
several known techniques. For example, the commonly 
used technique is to have the portable monitor the 
30 channel and determine which of several base stations 
it is closest to. It can then order the nearest BS to 
communicate with it. Another technique i6 to use a 
master control which receives information about the 
strength of the signal on the channel used by the 
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portable and controls the BS accordingly. Such 
techniques in themselves are known and do not form 
part of the invention. 

Such a protocol has several advantages. For 
5 instance, the location of the portable can be 
determined with high accuracy based on the received vc 
at the seven adjacent BS. Locating the portable can 
assist in the BS hand-off. A BS hand-off and a 
portable hand-off do not necessarily occur 
10 simultaneously, contrary to other prior art systems. 

In the present invention, when a portable roams from 
one cell X to an adjacent cell Y, a new vc is not 
required immediately. What is required is a BS hand- 
off, meaning that BS Y (associated with cell Y) must 
15 initiate transmission to the portable over the same 
vc, while the BS X (associated with cell X) must 
terminate its transmission to the portable. 

A BS hand-off occurs without the knowledge 
of the portable and can occur several times before a 
20 portable hand-off is required. A portable hand-off is 
required only when the CIR is below a certain level. 
In this case, the Mobile Telephone Switching Office 
(not shown) calls for a portable hand-off in 
accordance with known procedures. Reducing the 
25 portable hand-off rate reduces the probability of 
dropped calls . This is because a dropped call occurs 
either because the portable hand-off is not successful 
or because there are no available channels in cell Y. 
The present invention allows the use of post- 
30 detection diversity at the BS, and the use of dynamic 
channel allocation (DCA). 
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Dynamic Channel Allocation 

DCA is made possible by each BS having 
capability to transmit/receive more than the number of 
FDvc allocated to its cell, namely seven times the 
number of FDvc for a 5.86 MHz band and up to twenty- 
one times the number of FDvc for a 100 MHz as well as 
a 50 MHz band. The DCA protocol simply consists of 
borrowing as many FDvc as needed from the adjacent 
cells, up to a certain limit. The limit for the case 
when we employ a 6.18 Kbps vocoder / a 5.86 MHz band 
and 18 frames per cell is obtained as follows. The 
cell reuse pattern consists of 7 cells. Each cell is 
preassigned 144 FDvc. Assuming that at peak hours, 75 
FDvc are used on the average and 5 FDvc are reserved 
at all times, then we are left with 64 idle channels 
which represent the limit on the number of FDvc one 
can borrow from the cell. 

One should distinguish between the limit on 
the channels borrowed and the limit on the 
nonpreassigned channels a BS can use. For instance, if 
a call originates in- cell X and the portable roams 
into an adjacent cell Y where no preassigned cells are 
available, BS Y does not need to borrow immediately a 
new channel from an adjacent cell. It can use the 
original channel as long as the level of CIR is 
acceptable. If on the other hand, a portable wants to 
initiate a call in cell Y where all preassigned 
channels are used, BS Y can borrow a channel from an 
adjacent cell up to a limit of 64 channels per cell. 

The main advantage of DCA over Fixed Channel 
Allocation (FCA) is the increase in traffic handling 
capability. For FCA, a 7 cell pattern each with a 
preassigned 144 Fdvc can carry a total traffic of 
880.81 Erlang at 0.01 Blocking Probability (BP). For 
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DCA, a 7 cell pattern consists of 6 cells each with 80 
FDvc that can carry a total traffic of 392.17 Erlang, 
combined with one cell with 528 FDvc that can carry 
501.74 Erlang. The total traffic is therefore 893,91 
5 Erlang. This increase appears to be marginal (1.5%). 
However, if 501.74 Erlang are actually offered to one 
cell in the FCA system (with 14 FDvc/cell), while the 
six other cells carry 392.17/6 = 65.36 Erlang per 
cell, the BP at that busy cell 0.714 while it is 

10 negligible at the six other cells. The total blocked 
traffic (i.e. lost traffic) in the FCA system is then 
equal to (6x65.36x0.0 + 1x0.714x501.24) = 358.24 
Erlang. This represents a 0.4 average BP. If the DCA 
is allowed such a loss, its traffic handling capacity 

15 would increase to 1768.04 Erlang which represents a 
100% increase in traffic handling capacity over the 
FCA system, or equivalently a 160% increase in the 
number of available FDvc. The DCA system thus 
represents a marked improvement over the FCA system. 

20 

Voice Activation 

Voice activation is controlled by the BS 
according to techniques known in the art. At any 
instant during a conversation between a BS and a 
25 portable, there are four possibilities: 

1. BS talks while the portable listens. 

2. BS listens while the portable talks. 

3. BS and portable talk simultaneously. 

4. BS and portable listen simultaneously. 

30 The BS controls the voice activation 

procedure by allocating in cases 1, 3 and 4 three 
slots (frames 1.1, 1.2 and 1.3) to the BS and one slot 
the portable (frame 1) every four slots as shown in 
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Figure 9a. Likewise up to 21 portables may communicate 
with the base station in like fashion. 

In case 2, on receiving a signal from the 
portable, the BS allocates three slots (frames 1.1/ 
5 1.2 and 1.3} to the portable and one slot (frame 1) to 
the BS every four slots as shown in Figure 9b. 
Likewise, up to 21 other portables may communicate 
with the base station in like fashion. Consequently, 
instead of transmitting two full duplex voice frames 
10 over four slots as in Figure 4, voice activation 
allows us to transmit three full duplex voice frames 
over four slots. Hence, voice activation provides a 
50% increase in the number of available FDvc at the 
expense of increasing DSP complexity. 

15 

Capacity 

The capacity of Code Division Multiple 
Access (CDMA) may be defined as the number of half 
duplex voice channels (HDvc) effectively available 

20 over a 1.25 MHz band per cell. Based on such a 
definition, Table IV displays the capacity of analog 
FM and of the present system with a 6.18 Kbps vocoder, 
5.86 MHz band, 1 frame per cell and DCA. As shown in 
Table IV, the capacity of analog FM is 6 

25 HDvc/1 . 25MHz/cell while for the present system it is 
150 HDvc/1.25MHz/cell. The 6.25 MHz band consists of 
5.86 MHz plus two tail slots. When voice activation is 
used, the capacity of the present system is increased 
by 1,5 times to 225 HDvc/1 .25MHz/cell, a 38 fold 

30 increase over analog FM. 
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Local Area Networks 

The invention may also be applied to produce 
a 48 Mbps wireless LAN, which also satisfies the 
technical requirements for spread Bpectrum. 
5 For wireless LAN, wideband differential 

orthogonal frequency division multiplexing is again 
employed. The LAN will incorporate a plurality of 
transceivers, all more or less equal in terms of 
processing complexity, and possibly with identical 

10 components except for addresses. 

To implement wideband modulation for a LAN, 
a 26 MHz band is divided into 128 points, as shown in 
Figure 10, plus two tail slots of 1.48 MHz each within 
the 26 MHz band. Adjacent points are separated by 180 

15 KHz and each point, as with the application described 
above for a portable-base station, represents a D8PSK 
symbol. The transmitter components will be the same as 
shown in Figure 5b, with suitable modifications as 
described in the following, and will include an 

20 encoder. The output bits from the encoder are mapped 
onto the D8PSK symbols. 

The frame duration for the symbols is 
illustrated in Figure 11. A rectangular time domain 
window corresponding to a RC frequency domain pulse 

25 has a 5.55 jl/s duration, and includes a 25% roll-off 
and excess frame duration of 0.26 ys, making a total 
7.2 ys duration for the frame. 

For such a wireless local area network 
(LAN) , in which the transceivers are equal, the Time 

30 Division Duplex protocol is as illustrated in Figure 
12 (assuming there are at least a pair of 
transceivers) : 

1. A first transceiver transmits a signal (frame 0) 
over the entire frame. 
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2. A second transceiver receives the signal from the 
first transceiver and processes (analyzes) it. 

3. Based on the received signal, the second 
transceiver predistorts and transmits nine frames 

5 (frames 1 - 9 ) to the first transceiver immediately. 

Each transceiver has transmitter components 
similar to those illustrated in Figure 5b, with 
suitable modifications to the internal structure to 
allow the use of the particular frequency band and 

10 frame duration employed. 

The transmitter/receiver functional and 
structural block diagrams are shown in Figures 13a, 
13b and 13c for the exchange of data. Data is provided 
to an encoder 810 where the speech is digitized and 

15 coded to create bits of information. The bits are 
provided to the modulator 812 which turns them into 
D8PSK symbols, with three bit6 per symbol. The D8PSK 
symbols are then processed in the processor 814 which 
is described in more detail in Figure 14a. The output 

20 from the processor is then filtered in low pass filter 
816, upconverted to RF frequencies using local 
oscillator 818 and transmitted by antenna 820. 

In Figure 13b, the received Bignal at the 
base station is filtered in a bandpass filter 822, and 

25 down converted by mixing with the output of a local 
oscillator 824. The average power of the downcoverted 
signal is monitored by an initial power control 825 
that adjusts the average power to the specifications 
required by the sampler 826. The adjusted 

30 downconverted signal is then sampled in sampler 826 to 
produce bits of information. The bits are then 
processed in the deprocessor 828, described in more 
detail in Figure 14b. An estimate of the phaBe 
differential is taken in the channel estimator 830, as 
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described in more detail in relation to Figure 7 
above, and the estimated phase differential is 
supplied to a decoder/demodulator 832 to correct the 
received bits. The estimated phase differential is 
5 also supplied to a pre-distorter 834 in the 
transmitter. At the transmitter in the Base Station, 
the same blocks are incorporated as in the portable 
transmitter except that a pre-distorter is used to 
alter the phase of the D8PSK symbols to make the 

10 channel appear Gaussian (ideal) as opposed to a fading 
channel. The initial power control 825 also sends a 
signal to the pre-distorter 834 to adjust the 
transmitted power to an appropriate signal level for 
the sampler 826 in the first transceiver. It will be 

15 appreciated that a pre-distorter will be included in 
the first transceiver's transmitter but that it will 
not be operable, except when the first transceiver is 
operating as a base station. 

Figure 13c shows the functional blocks of 

20 the receiver of the first transceiver, which is the 
same as the receiver in the second transceiver except 
it does not include an estimator. The processor is 
illustrated in Figure 14a and 14c and the deprocessor 
in Figure 14b and 14c. The processor first inverse 

25 Fourier transforms the 128 D8PSK symbols output from 
the modulator. The transformed symbols are then 
triplicated as a group so that the total number of 
samples is tripled (see the left side of Figure 4c), 
with three consecutive groups each consisting of the 

30 128 transformed symbols. Next, the three groups are 
windowed by a Raised Cosine window with a roll-off of 
0.25 centered in the middle of the three groups. In 
other words, the processor takes D8PSK symbols in, 
pulse shapes them and inverse Fourier transforms them. 
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On the other hand, the deprocessor undoes what the 
processor did, i.e. it removes the pulse shaping, then 
Fourier transforms the received signal to obtain the 
original D8PSK symbols. The first two blocks in Figure 
5 14b are similar to the second two blocks in Figure 14a 

except for two differences as follows- In the first 
block shown in Figure 14b, the repeated groups of 
symbols are partially overlapped, as shown in Figure 
14c. In the second block, a rectangular window is used 

10 instead of the Raised Cosine to produce 128 output 
samples corresponding to the 416 input samples. 

The phase estimator is the same aB that 
shown in Figure 7, except that there are only 128 
input samples, and the same description applies. 

15 For both the LAN and cellular networks, the 

present system is designed to operate as a spread 
spectrum system preferably over such bands as are 
permitted, which at present are the 902 - 928 MHz 
band, 2.4 - 2.4835 GHz and 5.725 - 5.85 MHz - The 

20 carrier frequency in the local oscillator shown in 
Figures 5a, b and c will then be 915 MHz in the case 
of the 902 - 928 MHz band, and the frequencies used 
for modulation will be centered on this carrier 
frequency. 

25 

Alternative Embodiments 

A person skilled in the art could make 
immaterial modifications to the invention described 
and claimed in this patent without departing from the 
30 essence of the invention. 

For example, a system may consist of one or 
more central controllers (comparable to the Base 
Stations in the exemplary cellular system described) 
and some slave units (comparable to the portables). 
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39 

The slave unit executes the commands it receives from 
the central controller. The commands may be requesting 
the slave unit to transmit a receive acknowledge, a 
status code or information that the slave has access 
to. The command may also be to relay the command or 
the information to another slave unit. 
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BS COmnlftYitv of nrnrocenr/ 

de-processor/estimator 


vocoder 
rate 


complexity of processor/ 
de-processor for portable 


630 Mips 


18.77 Kbps 


1.35 Mips 


9.15 Kbps 


1.29 Mips 


6.18 Kbps 


1.26 Mips 



t Mips = Mega Instruction Per Second. 

1 1 1nstruction = 1 complex add.+1 complex multi.+1 storage. 



Table II 



Vocoder 
rate 


received 
SNR 


BER 


18.77 Kbps 


27.8 dB 


4X10* 3 


9.15 Kbps 


30.8 dB 


2X10" 3 


6.18 Kbps 


31.33dB 


2X10 3 



Table III 



2064975 



100 MHz, 1 frame/cell 


#of 
cells 


#of 
layers 


urban 

area, 

inKnrr 


indoor 
area - 
in Km 2 


CIR 

in dB 


126 


7 


24.6 


0.355 


37.9 



(a) 



50 MHz, 2 frames/cell 


#of 
cells 


#of 
layers 


urban 
area 9 
in Km 2 


indoor 
area- 
ta Km 2 


CIR 
in dB 


63 


5 


11.94 


0.172 


25.8 



(b) 



\ 



5.86 MHz, 1 8 frames/cell 


#of 
cells 


#of 
layers 


urban 
area, 
in Km 2 


indoor 
area, 
in Km 2 


CIR 
in dB 


7 


2 


1.43 


0.02 


13.8 



(c) 
Table IV 
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ccccccccccccccccccccccccccccccccccccccccccccccccc 



ccc 

c c 

c c 

C A LISTING FOR THE -PROCESSOR" C 

C WRITTEN IN FORTRAN C 

C C 



ccccccccccccccccccccccccccccccccccccccccccccccccc 
ccc 



Chh-^+++++++++++++++++hh-+++++++hh^++^+++^++++++++++++ 

C NAME : BOXNRM(ISTRM) - FUNCTION + 

C + 

C FUNCTION : SIMULATES NORMALY DISTRIBUTED DEVIATES + 

C ( POLAR METHOD OF BOX & MULLER IS USED) + 

C REF : + 

C + 

C CALLING PROG /SUBPROG: + 

C + 

C CALLED SUBPROG : URAND + 

C + 

C INPUTS : ISTRM - PASSED FROM THE CALLING ROUTINE + 

C + 

C OUTPUTS : STD. NORMALY DISTRIBUTED RANDOM DEVIATES 
+ 

C + 

C PRINCIPAL VARIABLES : + 

C + 

C DATE : 26 APRIL, 1985. + 

C + 

C+ I I I I I I I I ++++^+^+++++4^+++^++++++++++++++++++++++++++++ 

c 

FUNCTION BOXNRM(ISTRM) 

C 

C DECLARATION : 
C 

REAL URAND, V1, V2, R, FAC, GSET, BOXNRM 
INTEGER ISTRM, ISET 
SAVE GSET 

C 

DATA ISET/O/ 

C 

C SHOULD WE GENERATE A PAIR OF DEVIATES ? 
C 

IF(ISET.EQ. 0) THEN 

C 

C YES, WE DONT HAVE AN EXTRA DEVIATE HANDY, SO PICK 
C TWO UNIFORM NUMBERS IN THE SQUARE EXTENDING FROM 



2064975 



APPENDIX A 

Copyright § 1991 University Technologies Inc. 



C -1 TO +1 IN EACH DIRECTION 

C 

1 V1 =2. *URAND(ISTRM)-1. 

V2 = 2."URAND(ISTRM)-1. 

C 

C SEE IF THEY ARE IN THE UNIT CIRCLE, IF THEY ARE NOT GO TO 1 
C 

R = V1"2 + V2"2 
IF(R.GE. 1.) GOT01 

C 

C NOW MAKE THE BOX-MULLER TRANSFORMATION TO GET TWO 

NORMAL DEVIATES 

C 

FAC = SQRT(-2. * ALOG(R)/R) 

C 

C RETURN ONE AND SAVE ONE FOR NEXT TIME 
C 

GSET - V1 * FAC 
BOXNRM = V2 * FAC 

C 

C SET FLAG 
C 

ISET-1 

C 

ELSE 

C 

C WE HAVE AN EXTRA DEVIATE HANDY, SO RETURN IT 
C 

BOXNRM = GSET 

C 

C UNSET THE FLAG 
C 

ISET = 0 

C 

ENDIF 

C 

RETURN 
END 
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subroutine demod16(delta,cc) 
complex cc 

common /pie1/piover4,piover8,piover1 6 

if(delta.ge.31.*piover16.or .delta.lt. 1.*piover16)co(+1,+1) 
if(delta.ge. 1.*piover16.and.deIta.lt. 3/piover16)cc=(+1,-1) 
if(delta.ge. 3.*piover16.and.delta.lt. 5.*piover16)cc=(-1,-1) 
if(delta.ge. 5.*piover1 6. and.delta.lt. 7.*piover16)cc=(-1,+1) 

if(delta.ge. 7.*piover16.and.delta.lt. 9.*piover16)cc=(-1,+3) 
if(delta.ge. 9.*piover1 6.and.delta.It.1 1 /pioverl 6)cc=(-1 ,-3) 
if(delta.ge.1 1 .*piover16.and.delta.lt.13.*piover16)cc=(+1 1 -3) 
if (delta.ge.1 3. *piover1 6.and.delta.lt.1 5.*piover1 6)cc=(+1 ,+3) 

if(delta.ge.15.*piover16.and.delta.lt.17.*piover1 6)cc=(+3,+3) 
if(delta.ge.17.*piover16.and.delta.lt.19.*piover16)cc={+3,-3) 
if (delta.ge.1 9.*piover1 6.and.delta.lt.21 /piovert 6)cc=(-3,-3) 
if (delta.ge.21 .*piover1 6.and.delta.lt.23.*piover1 6)a>(-3,+3) 

if(deita.ge.23.*piover1 6.and.de)ta.lt.25.*piover1 6)cc=(-3,+1 ) 
if(delta.ge.25.*piover16.and.delta.lt.27/piov©r16)cc={-k3,+1) 
if(de!ta.ge.27.*piovert 6.and.delta.lt.29.*piover1 6)cc=(+3,-1 ) 
if(de!ta.ge.29.*piover1 6.and.delta.lt.31 ,*piover1 6)cc=(-3,-1 ) 

return 
end 
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subroutine demod4(delta,cc) 
complex cc 

common /piel /piover4,piover8,piovert 6 

if(delta.ge.0.*piover4.and.delta.lt.2.*piover4)cc=cmplx(+1 ,+1 ) 
if(delta.ge.2.*piover4.and.delta.It,4.*piover4)cc=cmplx(+1 ,-1) 
if(delta.ge.4.*piover4.and.delta.lt.6.*piover4)cc=cmplx(-1 ,-1 ) 
if(delta.ge.6.*piover4.and.delta.lt.8/piover4)cc=cmplx(-1 t +1) 

return 
end 
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subroutine demod8(delta,cc) 
complex cc 

common /piel/piover^pioverS.pioverl 6 

if(delta.ge.15.*piover8.or .delta.H 1.*piover8)co=(+1,+1) 
if(delta.ge. 1.*piover8.and,delta.lt. 3.*piover8)cc=(+1,-1) 
if(delta.ge. 3.*piover8.and.delta.lt. 5.*piover8)cc=(-1,-1) 
if(delta.ge. 5.*piover8.and.delta.lt. 7.*piover8)cc=(-1,+1) 

if(delta.ge. 7.*piover8.and.delta.lt. 9.*piover8)cc=(-3,+1) 
if(delta.ge. 9.*piover8.and.delta.It.1 1 .*piover8)cc=(+3,+1) 
if(de!ta.ge.11.*piover8.and.delta.lt.13.*piover8)cc=(+3,-1) 
if(delta.ge.13.*piover8.and.deita.it.15.*piover8)cc=(-3 l -1) 

return 
end 
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c SUBROUTINE FFT() - REFER TO THE BOOK BY 'GONZALEZ & WINTZ 
c 

c f is the input array to the fft 

c f is also the output array from the fft with +ve freq. at the begining 
c and -ve freq. at the end. 

c f 1 is the output array from the fft with -ve freq. at the begining 
c and +ve freq. at the end. 

c ind is the index array for f 1 . 

c if ix=0 the F is the input array, 
c if ix=1 the F1 is the input array. 



subroutine fft (ix.fjnd.fl, In) 

integer in, n, nv2, nm1 ( j, i, k, fe, Ie1, ip t ind(2048) 

real pi 

complex f(2048), f 1(2048) f u, w, t 

pU4.0*atan(1.0) 
n = 2**In 
nv2 = n/2 
nm1 = n-1 

if(ix.eq.1 )call deswap(f t f 1 ,n) 
j = 1 

do3i» 1, nm1 
if (i.ge.j) goto 1 
t-t(D 
f(i) = *(i) 
f(i) = t 
k= nv2 

if (k.ge.j) go to 3 



go to 2 
j«j+k 
do5l= 1, In 

le = 2**l 

Ie1 = le/2 

u = (1.0 t 0.0) 

w = cmplx( cos(pi/ie1 ), -sin(pi/le1 ) ) 
do5J=1,le1 
do4i==j,n ( le 
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ip = i + Ie1 
t-f(ip)*u 
f(ip) = f(i)-t 

4 f(i) = f(i) + t 

5 u = u * w 
do6i=1,n 

f(i) = f(i)/ float (n) 

6 ind(i) = i-n/2-1 

call swap(f,f1 ,n) 

return 

end 
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C""** * * " 

C Inverse fft transform 

C********** 

c 

c if ix=0 t f is the input array to the Ifft with +ve freq. at the begining 
c and -ve freq. at the end. 

c if ix=1 , f1 is the input array to the Ifft with -ve freq. at the begining 
c and +ve freq. at the end. 

subroutine tdft (ix ( xft,ind # xft1,nft) 

integer nft, i, mft,ind(2048) 

complex xft(nft),xft1(2048), buff(2048), buffi (2048) 

c 

mft - nint(alog10(float(nft))/alog1 0(2.0)) 
print \ ' nft and mft are:', nft, mft 

if(ix.eq.1)cali deswap(xft,xft1 ,nft) 

do 7 i= 1 , nft 

7 buff(i) = conjg(xft(i)) 

c fft() does 1 -D Fourier transform 

call fft(0 ( buff t ind t buff1 ( mft) 

do 8 i= 1 , nft 
xft(i) = nft*conjg(buff(i)) 

8 xft1(i) = nft*conjg(buff1(i)) 

return 
end 
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subroutine swap(f,f1 t n) 
complex f(2048), f 1(2048) 

do i=-n/2,-1 
f1(i+n/2+1) = f(n+i+1) 
enddo 



do i=0,n/2-1 
f1(i+n/2+1) = f(i+1) 
enddo 

return 
end 



subroutine deswap(f,f1,n) 
complex f(2048), f 1(2048) 

do i=-n/2,-1 
f(n+l+1)»f1(l+n/2+1) 
enddo 

do i=0,n/2-1 
f(i+1) = fl(i+n/2+l) 
enddo 

return 
end 
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real pi JWOPI,arRe(0:512),arlm(0:51 2) ( norm,p, 
& curphase(0:51 2),oldph,ph,psi(0:51 2),RC(-2048:2048) 
integer ind(2048) f counter 
complex 7(0:512)^(0:51 2) l c(0:512) i d(0S12) l cc 
complex F(2048),F1(2048) 
character string*50 
common /pie/pi 

common /pie1/piover4 l piover8,piover1 6 
common /dec/no 

open(1 OO.file^'transmitted') 

open(200,file-'received*) 

open(400,file= , BER l ) 

open (800 ( f ile=envelope') 

open^lO.file^RC) 

open(350 ( file= , data , ) 

open(351 ,file= random') 

pi = 4. * atan(1.0) 

piover4 = pi/4. 

piover8 = pi/8. 

piover16 « pi/16. 

TWOPU2.*pi 

read(M02)string 

write(\102)string 

read(* *)number 

read(Y)M 

read(Y)MM 

read(Y)no 
102 format(a50) 
C open(1 01 ,fite=string) 

C numb=255 
LL= MM/2 -1 
counter = 0 

call raisedcos(MM,RC) 
C doi»1,numb 

c***"" READING * # 

C read(1 01 , # }p l arRe(i),arlm(i) 

C norm »( 900. /p)"2 

C arRe(i) = arRe(i)/norm 

C arlm(i) = arIm(i)/norm 

C z(i)=cmplx(arRe(i),arlm(i)) 

C write(800,*)l.sqrt(real(Z(i)*conjg(Z(i)))) 

C enddo 

C 2(0) =cmp!x(arRe(1 ),arlm(1 )) 

C z(MM)=cmplx(arRe(numb),arlm(numb)) 
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do ii=1, number 
C""" FIRST LOOP 

curphase(O) = urand(3)*TWOPI 

d(0) = cmplx(cos(curphase(0)),sin(curphase(0))) 

C* STARTING THE DATA TRANSMISSION 
C**"" DATA GENERATION •«••»••••—••••••••«• 

do10l=1,MM,1 
if(M.eq.4)call modu4(cc,delta) 
if(M.eq.8)call modu8(cc,delta) 
if(M.eq.16)call modu16(cc,delta) 
write(500,*)l,cc,delta 

c(l) = cc 

psi(l) = urand(4)*TWOPI 
write(35V)psi(l) 

curphase(l) = curphase(l-1 ) + delta + psi(l) 
d(l) = cmplx(cos(curphase(l)),sin(curphase(l))) 
enddo 

C""" MODULATOR 

do k=0,MM,1 
write(100,*)k,sqrt(real(d(k)*conjg(d(k}))) 
C r(k) = d(k)*z(k) 
r(k) = d(k) 
write(451,*)k,r(k) 
enddo 

C "*"* FIND THE TIME RESPONSE 
j=0 

dok=1,MM,1 
H+1 

Fl0) = r(k-1) 
enddo 

CALLtdft(1,F,ind,F1,MM) 

open(300,file='time_RCSF) 
open(301 ,file='timeO') 
openpoa.file-'timejvlOSF) 

call window(MM,LL,RC,F1 .counter) 

C goto101 

c ...... DEM ODULATOR "***" 

oldph = atan2(aimag(r(0)),real(r(0))) 



9 
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do k=1,MM-1 
ph = atan2(aimag(r(k)),real(r(k))) 
delta = ph - oldph - psi(k) 
if(delta.gt.TWOPI)deita=delta-TWOPI 
if(delta.gt.TWOPI)de!ta=delta-TWOPI 
if(delta.gt.TWOPI)delta=delta-TWOPI 
if(deltalt.0.)delta=delta+TWOPI 
if(delta.lt.0.)delta=delta+TWOPi 
if(deita.lt.0.)deita=delta+TWOPI 

if(M.eq.4)cail demod4(delta.cc) 
if(M.eq.8)call demod8(delta,cc) 
if (M.eq. 1 6)call demodl 6(delta,cc) 
write(501,*)k,cc,delta 

if(real(cc).ne.real(c(k)))errorRe=errorRe+1 
if(aimag(cc).ne.aimag(c(k)))errorlm=errorlm+1 

wrfte(400.*)k/ \real(cc),real(c(k)), 
& 1 \aimag(cc) f aimag(c(k)) 

oldph = ph 
enddo 

printYerrorRe^errorRe,* 7errorim=\errorlm 
10 continue 
101 continue 

enddo 

rewind(301) 

call quantize(LL,counter,prod) 



stop 
end 
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subroutine modu16(cc,delta) 
complex cc 

common /pie1/piover4,piover8 l piover1 6 

s1 - urand(1) 

s2 = urand(2) 

s3 = urand(1) 

s4 = urand(2) 

if(s1.ge.0.5)s1=+1 

if($1.lt.0.5)s1=-1 

if(s2.ge.0.5)s2=+1 

if(s2.lt.0.5)s2=-1 

if(s3.ge.0.5)s3=+1 

if(s3.lt.0.5)s3=-1 

if(s4.ge.0.5)s4=+1 

if(s4.lt.0.5)s4=-1 

cc = cmplx(s1+2*s3,s2+2*s4) 

write(350,*)cc 

if(cc.eq.(+1 1 +1 ))delta=0.*piover8 
if(cc.eq.(+l ,-1 ))delta=1 .*piover8 
if(cc.eq.(-1,-1))delta=2.*piover8 
if(cc,eq.(-1 ,+1 ))delta=3.*piover8 

if(cc.eq.(-1,+3))delta=4.*piover8 
if(cc.eq.(-1,-3))delta=5.*piover8 
if(cc.eq.(+1 ,-3))delta=6.*piover8 
if(cc.eq.(+1,+3))delta=7.*piover8 

if(cc.eq.(+3,+3))delta=8.*piover8 
if(cc.eq.(+3r3))deIta=9.*piover8 
ifjcc.eq.l-SrSJJdelta^lO/piovera 
if(cc,eq.(-3,+3))deIta=1 1 /pioverS 

if(cc.eq.(-3,+1))delta==12.*piover8 
if(cc.eq.(+3,+1 ))delta=1 3.*piover8 
if (cc.eq.(+3,-1 ))delta=1 4.*piover8 
if (cc.eq.(-3,-1 ))delta=1 5.*piover8 

return 
end 
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subroutine modu4(cc,delta) 
complex cc 

common /pie1/piover4,piover8,piovert 6 

s1 = urand(1) 
s2 = urand(2) 
if(s1.ge.0.5)s1=+l 
if(s1.lt.0.5)s1=-1 
if(s2.ge.0.5)s2=+1 
if(s2.lt.0.5)s2=-1 
cc = cmplx(s1 ,s2) 
write(350,*)cc 

if(cc.eq.(+1 ,+1 ))delta=1 .*piover4 
if(cc.eq.(+1 ,-1 ))delta=3/piover4 
if(cc.eq.(-1 ,-1 ))de!ta=5.*piover4 
if(cc.eq.(-1,+1))deIta=7.*piover4 

return 
end 
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subroutine modu8(cc,delta) 
complex cc 

common /pie1/piover4,piover8 ) piover1 6 

s1 = urand(1) 

s2 = urand(2) 

s3 = urand(1) 

if(s1.ge.0.5)s1=+1 

if(s1.lt.0.5)s1=-1 

if(s2.ge.0.5)s2=+1 

if(s2.lt.0.5)s2=-1 

if(s3.ge.0.5)s3=+1 

if(s3.lt.0.5)s3=-1 

cc » cmplx(s1+2*s3,s2) 

pp « urand(4)*TWOPI 

write(350,*)cc 

if(cc.eq.(+1 ,+1 ))delta=0.*piover4 
if(cc.eq.(+1 ,-1 ))delta=1 .*piover4 
if(cc.eq.(-1 ,-1 ))delta=2.*piover4 
if(cc.eq.(-1,+1))delta=3.*piover4 

if(cc.eq.(-3,+1))delta=4.*piover4 
if(cc.eq.(+3,+1))delta=:5.*piover4 
if (cc. eq. (+3,-1 ))delt a=6 , *piover4 
if(cc.eq.(-3,-1))delta^7.*piover4 

return 
end 
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subroutine quantize(LL ( counter,prod) 
integer Fre.Fim.counter 

xmax=-100. 
xmin=100. 
do j=1 .counter 

read(301,*)i,arRe,arim 

if(arRe.gt.xmax)xmax=arRe 

if(arlm.gt.xmax)xmax=arlm 

if(arRe.lt.xmin)xmin=arRe 

if(arlm.lt.xmin)xmin=arlm 
enddo 

if(abs(xmax).gt.abs(xmin))xxx=abs(xmax) 
if(abs(xmax).lt.abs(xmin))xxx=abs(xmin) 
prod = 128./xxx 

sumRe = 0. 
sumlm « 0. 
surRe = 0. 
surlm = 0. 
rewind(301) 
do j=1 .counter 
read(301,*)i,arRe.arlm 
Fre= 128 + prod*arRe 
Fim=s 128 + prod*arlm 
C if(Fre.gt.255)Fre=255 
C if(Fre.lt.0)Fre=0 
C if(Fim.gt.255)Fim=255 
C if(Fim.It.0)Fim=0 

write(401,*)i,Fre,Fim 
C write(402,*)j,Fre 
C write(403 ( *)j,Fim 

sumRe « sumRe + (arRe-(Fre-128.)/prod)**2 
sumlm « sumlm + (arlm-(Fim-1 28.)/prod)**2 
surRe = surRe + (arRe)**2 
surlm = surlm + (arlm)**2 
enddo 

printVquantization noise=',(sumRe + sumlm)/(surRe + surlm) 
print*,prod 

return 
end 
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subroutine raisedcos(MM.RC) 
real RC(-2048:2048).pi 
common /pie/pi 

do i=-MM/2-MM/8,-MM/2+MM/8 

t = -1loat(i)/float(MM) 

RC(i)=(0.5*(1.-sin(4.*pi*(t-0.5)))) 

write(310.*)i,RC(i) 
enddo 

do i=-MM/2+MM/8+1,MM/2-MM/8-1 
t = float(i)/float(MM) 
RC(i)=1 . 

write(310,*)i,RC(i) 
enddo 

do i=MM/2-MM/8,MM/2+MM/8 

t « float(i)/lloat(MM) 

RC(i)=(0.5*(1.-sin(4.*pi*(t-0.5)))) 

write(310,*)i,RC(i) 
enddo 

return 
end 
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subroutine demod4(delta,wRe,wlm) 

y = cmplx(cos(delta),sin(delta)) 
if(real(y).ge.O)then 

wRe=+1. 
else 

wRe=-1. 
endif 

if(aimag(y).ge.O)then 

wlm=+1. 
else 

wlm=-1. 
endif 

return 
end 
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C^+++^^++++++++^+++++++++++++++++-H-+++++++++++++++++++++ 

C NAME : URAND(ISTRM) • FUNCTION + 

C + 

C FUNCTION : GENERATES STANDARD UNIFORMLY DISTRIBUTED 

+ 

C RANDOM NUMBERS . + 

C ( USES THE RECURSION: SEED(ISTRM)= + 

C 1 6807 * SEED(ISTRM) MOD (2"(31 ) - 1 ) ) + 

C SOME COMPILERS REQUIRE THE DECLARATION: 

+ 

C INTEGERS ISTRM. K1 + 

C REF : + 

C + 

C CALLING PROG /SUBPROG: + 

C + 

C CALLED SUBPROG :NIL + 

C + 

C INPUTS : ISTRM -THE STREAM NUMBER + 

C 0 < SEED(ISTRM) < 21 47483647 + 

C + 

C OUTPUTS : STD. UNIFORMLY DISTRIBUTED RANDOM NUMBERS 

+ 

C NEW VALUE OF SEED(ISTRM) + 

C + 

C PRINCIPAL VARIABLES : + 

C + 

C DATE : 26 APRIL, 1985. + 

C + 

C4HHHH^^++-HHH-+4H^++++++++4+++^ 

c 

FUNCTION URAND(ISTRM) 

C 

C DECLARATION : 
C 

INTEGER SIZE3 
PARAMETER(SI2E3 = 20) 
REAL URAND 

INTEGER SEED(SIZE3).ISTRM,ISTRG,K1 

C 

save seed 

data seed(1)/456789/ 
data seed(2)/1 037625857/ 
data seed(3)/i 02899405/ 
data seed(4)/535417/ 
data seed(5)/2H 229747/ 
data seed(6)/2952853/ 
data seed(7)/2987570/ 
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data seed(8)/45376256/ 
data seed(9)/1 7890857/ 

C 

C GENERATE A U(0.1 ) VARIATE 
C 

ISTRG = ISTRM 

K1 « SEED(ISTRG)/1 27773 

SEED(ISTRG) = 16807 * ( SEED(ISTRG) - K1 * 127773) - K1 * 2836 
IF( SEED(ISTRG) .LT. 0) SEED (ISTRG) = SEED(ISTRG) + 2147483647 
URAND = FLOAT(SEED(ISTRG)) * 4.65661 2875E-10 

C 

RETURN 
END 
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subroutine window(MM,LL ( RC,F,counter) 

complex F(2048) l F0(-2048:2048) f Fm(-2048:2048),Fp(-2048:2048), 

Fsum(-2048:2048),FRC(-2048:2048) 
real RC(-2048:2048) 
integer counter,count1 
common /dec/no 
save countl 
data countl /O/ 

countl = countl + 1 
do i=-MM-LL,MM+LL+1 

FO(i)-(a.o.) 
Fm(i) = (0..0.) 
Fp(i) = (0..0.) 
enddo 

do i=-LL,LL+1 

F0(i) =F(i+LL+1) 

Fm(i-MM) = F0(i) 

Fp(i+MM) - F0(i) 

write(304, # )i,F0(i) 
enddo 

do i=-MM-LL,MM+LL+1 
Fsum(i) = F0(i) + Fm(i) + Fp(i) 

wiite(302,*)i,sqrt((real(Fsum(i)))"2+(aimag(Fsum{i)))**2) 
enddo 

do i=-LL - (LL+1)/4-8,-LL - (LL+1J/4-1 

counter = counter + 1 

FRC(i) - cmplx(0. ( 0.) 

write(301 ( *)counter,0.,0. 
enddo 

do i=-LL- (LL+1)/4,LL+1 + (LL+1)/4 
counter = counter + 1 
FRC(i) = Fsum(i)*RC(i) 

write(300,*)counter ( sqrt((real(FRC(i)))**2+(aimag(FRC(i)))**2) 
write(301/)counter,real(FRC(i)) t aimag(FRC(i)) 
write(202/)counter l real(FRC(i)j 
writef^OS.^counter.aimagtFRCfl)) 
enddo 

do i=LL+1 + (LL+1)/4+1,LL+1 + (LL+1)/4+8 
counter = counter + 1 

FRC(i) = cmplx(0.,0.) 
write(301 ,*)counter,0. f 0. 
enddo 
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if(count1.eq.no)then 

do i=-LL- (LL+1)/4-8,LL+1 + (LL+1)/4+8 

counter = counter + 1 

FRC(i) m cmplx(0.,0.) 

write(301 ,*)counter,0..0. 

write(202,*)counter,rea!(FRC(i)) 

write(203,*)counter,aimag(FRC(i)) 
enddo 
endif 

return 
end 
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FFLAGS - 

FFILES = main.f urand.f fft.f window/f raisedcos.f quantize.! \ 
modu4.f moud8.f modu16.f demod4.f demod8.f demod16.f 

OFILES « main.o urand.o fft.o window.o raisedcos.o quantize.o \ 
modu4.o modu8.o modu16.o demod4.o demod8.o demod16.o 

a.out: ${OFILES} 

f77 ${FFLAGS} ${OFILES} 
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../../.7../../HF_LF/TWOCHANNEL/900/test.1 
9 

16 
64 
9 
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ccaocnjjxxxxxixxxxxjjxjoooca 

c c 

C WRITTEN IN FORTRAN C 

c c 

C A LISTING FOR THE "DEPROCESSOR" C 
C C 

cccccccccccccccccccccccccccccccccccccccccccccccc 
cccc 



subroutine demod16(delta,cc) 
complex cc 

common /piel /piover4,piover8,piover1 6 

if(delta.ge.31.*piover16.or .delta.lt. 1 /plover! 6)ce={+1,+1) 
if(delta.ge. 1.*piover16.and.delta.lt. 3.*piover16)cc=(+1.-1) 
lf(delta.ge. 3.*piover16.and.deita.lt. 5.*piover16)cc=(-1,-!) 
if (delta.ge. 5.*piover1 6.and.delta.lt. 7.*piover16)cc=(-1,+1) 

if(delta.ge. 7.*piover16.and.delta.lt. 9.*piover16)cc=(-1 ( +3) 
if(delta.ge. 9.*piover1 6.and.delta.lt.1 1 /piover! 6)cc=(-1 ,-3) 
if (delta.ge. 1 1 .*piover1 6.and.deita.lt.1 3/pioverl 6)cc=(+! ,-3) 
if(delta.ge.1 3.*piover1 6.and.delta.lt.l 5.*piover1 6)cc=(+1 ( +3) 

if (delta.ge. 1 5.*piover1 6.and.delta.lt. 1 7.*piover1 6)cc=(+3,+3) * 
if (delta,ge. 1 7.*piover1 6.and.delta.lt. 1 9.*piover1 6)cc=(+3,-3) 
if(delta.ge.1 9.*piover1 6.and.delta.lt.21 .*piover1 6)cc=(-3,-3) 
if(delta.ge.21 /pioverl 6.and.delta.lt.23.*piover1 6)cc=(-3,+3) 

if(de!ta.ge.23. # piovert 6.and.delta.lt.25.*piover1 6)co=(-3,+1 ) 
if(delta.ge.25.*piover1 6.and.delta.lt.27.*piover1 6)co(+3,+1) 
if(delta.ge.27.*piover1 6.and.delta.lt.29.*piover1 6)cc=(+3,-1 ) 
if (delta.ge.29.*piover1 6.and.delta.lt.31 .'piovert 6)cc=(-3.-1 ) 



return 
end 
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subroutine demod4(delta,cc) 
complex cc 

common /pie1/piover4,piover8,piover16 

if(delta.ge.0.*piover4.and.delta.lt.2.*piover4)cc=cmplx(+1 ,+1 ) 
if(delta.ge.2.*piover4.and.delta.lt.4.*piover4)cc=cmplx(+1 ,-1 ) 

if(delta.ge.4.*piover4.and.delta.lt.6.*piover4)cc=cmplx(-1 ,-1 ) 
if(delta.ge.6.*piover4.and.delta.lt.8.*piover4)cc=cmpix(-1 ,+1 ) 

return 
end 
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subroutine demod8(delta,cc) 
complex cc 

common /pie1/piover4,piover8,piover16 

if(de!ta.ge.15.*piover8.or .delta.lt. 1/piover8)cc=={+1,+1) 
if(delta.ge. 1.*piover8.and.delta.lt. 3.*piover8)cc=(+1,-t) 
if(detta.ge. 3.*piover8.and.delta.lt. 5.*piover8)c&=(-1,-1) 
if(delta.ge. 5.*piover8.and.delta.lt. 7.*piover8)cc=(-1,+1) 

if(deita.ge. 7.*piover8.and.delta.lt. 9.*piover8)cc=(-3,+1) 
if(delta.ge. 9.*piover8.and.delta.lt.1 1 .*piover8)cc=(+3,+1) 
if(d-9lta.ge.1 1 .*piover8.and.delta.lt.1 3.*piover8)c<H+3.- 1 ) 
if (delta.ge.1 3.*piover8.and.delta.lt.1 5/piover8)cc=(-3,-1 ) 

return 
end 
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subroutine dewindow{MM,numb,F) 
integer shift 

complex F(2048).FO(-2048:2048),Fm(-2048:2048), 
& Fsum(-2048:2048),Fp(-2048:2048) ( fprime 
common /SHHH/shift 
common /SDHH/prod 

do M.numb 

readCIOVJj.arRe.arlm 

F(i) « cmplx(arRe-shift,arlm-shift) 
enddo 

if(coun1.eq.no)then 
do i=1 ( numb 

read(1 01 ,*)j t arRe,arlm 
enddo 
endif 

LUMM/2-1 

do i=-MM-numb/2+1,MM+numb/2 

F0(i)-(0.,0.) 
Fm(i) = (0..0.) 
Fp(i) = (0.,0.) 
enddo 

do i=-numb/2+1,numb/2 

FO(i) -F(i+numb/2) 

Fm(i-MM) = F0(i) 

Fp(i+MM) = F0(i) 
enddo 

do i=-MM-numb/2+1 f MM+numb/2 

Fsum(i) m F0(i) + Fm(i) + Fp(i) 
enddo 

sum = 0. 
sup = 0. 
do i=-LL,LL+1 

write(303 ( *)i ( sqrt{(real(Fsum(i)))"2 
& +(almag(Fsum(i))r2) 

F(i+LL+1) = Fsum(i) 

write(304,*)i,Fsum(i) 

C read(604,*)j,Fprime 

C sum = sum + conjg(Fsum(i)/prod-Fprime)* 

C & (Fsum(i)/prod-Fprime) 

C write(602,*)i l abs(real(Fsum(i)/prod)) 
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C & -abs(reaI(Fprime)) 

C write(603,*)i 1 abs(aimag(Fsum(i)/prod)) 

C & -abs(aimag(Fprime)) 

C write(605, # )i,Fsum(i)/prod-Fprime 

C 

C sup = sup + conjg(Fprime)*(Fprime) 
enddo 

C printYquantization error= , l sum,sup ( sum/sup 

return 
end 
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c SUBROUTINE FFT() - REFER TO THE BOOK BY 'GONZALEZ & WINTZ 
C 

c f is the input array to the fft 

c f is also the output array from the fft with +ve freq. at the begining 
c and -ve freq. at the end. 

c f 1 is the output array from the fft with -ve freq. at the begining 
c and +ve freq. at the end. 

c ind is the index array for f 1 . 

c if ix=0 the F is the input array, 
c if ix=1 the F1 is the input array. 



subroutine fft (ix,f,ind,f1 , In) 

integer In, n, nv2 ( nm1, j, i, k, le, ie1, ip,ind(2048) 

real pi 

complex f(2048), f1(2048),u ( w, t 

c 

pi = 4.0 * atan(1 .0) 
n = 2**ln 
nv2 = n/2 
nm1 = n-1 

if(ix.eq.1)cali deswap(f,f1 t n) 
J-1 

do 3 i« 1, nm1 
if (i.ge.j) go to 1 
t = f(j) 

f(i) = t 

1 k = nv2 

2 if (lege j) go to 3 



go to 2 

3 j = j+k 

do 5 1=1, In 
le m 2**l 
Ie1 = le/2 
u»'(1.0 f 0.0) 

w = cmplx( cos(pi/le1), -sin(pi/le1) ) 
do5j=1,lel 
do 4 1= j, n, le 
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ip = i + Ie1 
t-f(ip)*u 
f(ip) = f(i)-t 

4 1(i) = f(i) + t 

5 u = u * w 
do6i=1,n 

f(i) = f(i)/ float (n) 

6 ind(i) = i-n/2-1 

call swap(f,f1,n) 

return 

end 
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C ****** **** 

C Inverse fft transform 

c 

c if ix=0, f is the input array to the Ifft with +ve freq. at the begining 
c and -ve freq. at the end. 

c if ix=1 f f1 is the input array to the Ifft with -ve freq. at the begining 
c and +ve freq. at the end. 

subroutine tdft (ix.xftjnd.xftl ,nft) 

integer nft. i, mft ( ind(2048) 

complex xft(nft) f xft1(2048), buff(2048), buffi (2048) 

c 

mft = nint(alog10(fioat(nft))/a;og1 0(2.0)) 
print *, ' nft and mft are:', nft, mft 

if(ix.eq.1)call deswap(xft,xft1 ,nft) 

do 7 i= 1 , nft 

7 buff(i) = conjg(xft(i)) 

c fft() does 1 -D Fourier transform 

call fft(0,buff t ind,buff1,mft) 

do 8 i= 1 , nft 
xft(i) - nft*conjg(buff(i)) 

8 xft1(i) = nft*conjg(buff1(i)) 

return 
end 
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subroutine swap(f,f1,n) 
complex f(2048), f1 (2048) 

do i=-n/2,-1 
f1(i+n/2+1) = f(n+i+1) 
enddo 

do i=0,n/2-1 
f1(i+n/2+1) = f(i+1) 
enddo 

return 
end 
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subroutine deswap(f ,f 1 , n) 
complex f(2048), f 1(2048) 

do i=-n/2,-1 
f(n+i+1) = fl (i+n/2+1) 
enddo 

do i=0,n/2-1 
f(i+1) = f1 (i+n/2+1) 
enddo 

return 
end 
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real piJWOPI,oldph,ph,psi(0:512) 

complex r(0:512),c(0:512),cc ( cl 

complex F(2048),F1{2048) 

integer ind(2048),shift 

character string*40,string1 *40,string2*40 

common /pie/pi 

common /pie1/piover4,piover8 ; piover1 6 
common /SHHH/shift 
common /SDHH/prod 

open(l OCfile^'transmitted') 
open(200 t file='received') 
open^O^AIe^'BER 1 ) 
openlsOO.file^envelope') 
openfSIO.file^RC 1 ) 
pi = 4. *atan(1.0) 
piover4 « pi/4. 
piover8 = pi/8. 
piover16 = pi/16. 
TWOPi » 2. * pi 
read(*,102)string 
write(M02)string 
read(M02)stringl 
write (M02)string1 
read(*,102)string2 
write(*,102)string2 
read(*/)number 
read(Y)M 
read(Y)MM 
read(Y)ni 
102 format(a40) 

open(101 f fiie=string) 
open(350,file=string1 ) 
open(351 ,ftle=string2) 
open^.file^./fort.SO^) 
open(704 ( file= , ../fort.501') 
shift = 0 
prod = 1. 

if(string.eq.\./foit.401')shift=1 28 

if (string.eq.\./fort.401 ')read(306,*)prod 



numb=MM*1.25+ 16 



c"***" READING 
sum = 0. 
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do ii=1, number 

call dewindow(MM,numb,F) 

call fft(0,F,ind.F1,ni) 
j-0 

dok=0,MM-1,1 
H+1 

r(k) = F10) 
write(451.*)k.r(k) 
read(350.*)c(k+1) 
read(351,*)psi(k+1) 
enddo 

C""" DEMODULATOR «——••••••• 

oldph = atan2(aimag(r(0)),real(r(0))) 

dok=1.MM-1 
ph = atan2(aimag(r(k)),real(r(k))) 
delta = ph - oldph • psi(k) + pi 
if(delta.gt.TWOPI)delta=delta-TWOPI 
if(delta.gt.TWOPI)delta=delta-TWOPI 
if (delta.gt.TWOPI)delta=delta-TWOP I 
if(delta.lt.O.)delta=delta+TWOPI 
if(delta.lt.0.)delta=delta+TWOPI 
if(delta.lt.0.)delta=delta+TWOPI 

if(M.eq.4)call demod4(delta,cc) 

if(M.eq.8)call demod8(delta,cc) 

if(M.eq.16)call demodl 6(delta,cc) 

write(50V)k,cc,delta 

read(704,*)l,cl,deltap 

diff = delta-deltap 

if(diff.gt.+pi)diff=diff-TWOPI 

if(diff.lt.-pi)diff=diff+TWOPI 

if(k.ne.MM/2.and.k.ne. MM/2+1 )sum = sum + diff**2 

write(606,*)k,diff 

if(real(cc).ne.real(c(k)))errorRe=errorRe+1 
if(aimag(cc).ne.aimag(c(k)))errorlm=errorlm+1 

write^Ou.*^. 1 , ,real(cc),real(c(k)). 
& ' ',aimag(cc),aimag(c(k)) 

oldph = ph 
enddo 

print*,'errorRe=',errorRe,' , ,'errorlm=',errorlm 
print*,'sum=',sum/(loat(MM-3) 
1 0 continue 
enddo 
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stop 
end 
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FFLAGS • 

FFILES = main.f fft.f dewindow.f demod4.f demod8.f demod16.f 

OFILES = main.o fft.o dewindow.o demod4.o demod8.o demod16.o 

a.out: ${OFILES} 

f77 ${FFLAGS} ${OFILES} 
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.Vfoit.401 
../data 
../random 
9 

16 
64 
6 
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SOFTWARE FOR SIMULATION ON DSP32C 

The following software program was developed and run using D3SIM of AT&T. It 
simulates the operation of the equalizer on a OSP32C. 

^include <madi.h> 
ttincludc <stdio.h> 
^include <libap.h> 

main program 

**+*»*******************************************»**********++***»***#**^ 

double twopi=6.283 185372; 

/* List the original transmitted symbols */ 

intjj2[514]={ -1 , -1, 

1 . 1. 

1 . -1. 

1 , -I. 

1 . -1. 
-1 . -1, 
-1 . 1. 

1.1. 
-1 , 1. 

1 . -1. 
-1 . -1. 

1 . 1. 
-I . -1. 
-1 . -1. 
-1 . 1. 
-1 . -1. 
-1 . 1. 
-1 . -1. 

1 . 1. 
-1 . -1. 
1 . -I. 
-1 . -1. 

1 . -1. 

1 . -1. 
-1 . 1. 
-1 . I. 

1 . -1. 

1 . -1. 
-1 . -1. 

1 . -1. 
•1 . -1. 
-1 . -1. 
-1 , 1. 
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-1, 
1, 
-1, 
-1, 

1, 

1, 

1. 
-1. 

1, 
-1, 

1. 
-1, 

1. 

1. 

1. 
-1. 
-1, 
-1. 

1, 
-1, 

1, 
-1. 

1, 

1, 
-1. 
-1. 

1. 
-1, 
-1, 

1. 
-1. 
-1, 

1. 
-1. 
-1, 
-1. 

1, 
-1. 
-1, 
-1. 
-1. 
-1. 
-1, 

1. 
-I. 
-1, 

1. 
-1. 
-1, 

1, 

1. 
-1, 
-1. 



2 



£064975 



APPENDIX C 

Copyright 8 1991 AGT Limited 



-1, 
-1, 
-1. 

1. 
-1, 
-1. 

1. 
-1, 
-1. 
-1. 
-1, 

1. 
•1. 

1. 

1. 

1. 
-1. 
-1. 
-1, 

-1, 

1, 

1. 

1. 

1. 

1. 

1. 

1. 
-1. 
-1, 
-1. 

1. 
-1, 
-1, 
-1, 

1, 
-1. 
-1. 
-I. 

1. 

1. 

1, 
-1. 

1, 

1, 
-1. 

1. 

1. 
-1, 

1, 
-1, 
-1. 

1. 
-1. 
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. 1, 
, -1, 
, -1. 
. -1, 
. -1, 
. 1, 
, -1. 
. 1, 
. 1. 
. -1. 
. -1, 
. 1, 
. 1, 
, -1, 
. 1. 
. -1, 
. -1, 
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/* List the received I and Q */ 
double z[514]={ 
1.8826647996902 ,0.57904273271561. 

-0.98167502880096 , -1.7151488065720. 

-1.9221851825714 . -0.45452463626862, 

-1.6370518207550 , 1.0898391008377. 

-0.31615784764290 , 1.9248703718185. 

-1.5235337018967 , 1.1806204319000. 
0.16620002686977 , -1.8899723291397, 

-1.3761224746704 , 1.2521643638611. 
7.1681058034301E-03 , -1.8177816867828. 

-1.3031834363937 .-1.1972546577454, 

-0.15820230543613 , -1.7096272706985. 

-0.99022841453552 , 1.3329787254333. 
-1.567842841 1484 , 0.32700437307358. 
-1.3414373397827 , -0.75913363695145, 
1.3957004547119 , -0.49621558189392. 
1.3290159702301 , 0.50875544548035. 
-1.1973139047623 , 0.66275292634964, 
-1.2967063188553 , -0.24445547163486. 
0.82353335618973 , 0.97751879692078, 
2.7966260910034E-02 . 1.2459809780121, 
0.97553628683090 , 0.74172759056091, 
0.30242365598679 , 1.1787307262421, 
1.1158645153046 , 0.49576562643051, 
1.0971860885620 . -0.57270908355713. 
-0.2456972S99029S , 1.2418062686920, 
-0.83266228437424 , -1.0038361549377, 
-1.3508954048157 , 2.3582908324897E-03. 
-0.90133684873581 , 1.0763428211212, 
-1.4409677982330 , 024240809679031, 
-0.79242151975632 , 1.2981982231140, 
-1 .5102118253708 , 0.4687696993351O, 
-1.4932200908661 , -0.67981052398682, 
1 .5572143793106 . -0.67623662948608, 
1.6570891141891 , 036612873077393, 
-13809999704361 , 0.86023199558258. 
-1.7863020896912 , -0.45382413268089, 
-1.5810590982437 , 1.0169448852539. 
0.34510406851768 , -1.8782730102539, 
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1 1434341669083 . 1.5573689937592, 

- 5104018449783 . 1.2377147674561. 
1 9451805353165 , 0.14569547772408. 
1 90SS004684448 . 1.4411246776581, 

-1 3509848117828 , 1.3267296552658, 

i5?7233552933 ±™gfflflg™ a > 
13225439786911 ,-1.2418869733810, 
ffloiHS -9.0550981462002E-02. 
1 1161555051804 . -1.2882466316223. 
1 6323974132538 . -0.1504679173231 1. 
097648859024048 .-1.2267235517502, 

-0.82589977979660 . \-}4^™g%' 
1 3004480600357 . -0.24329216778278. 
oSlSSSlMS .-1-0373107194901. 
0 27767503261566 . 1.1074353456497. 
0 50518971681595 -0.91855829954147 
M0428864955902 -0.30539751052856. 
079054325819016 . 0.34204551577568. 
-0 69724625349045 0.32767096161842. 
-0 18177013099194 . 0.65859514474869. 
0 34578496217728 0.49258133769035. 
0 52803617715836 . -2.7841679751873E-02. 
:0 29641786217690 . 0.36104345321655. 
-64Ko594472885 0.11642142385244. 
011455089598894 -0.37470200657845. 
0 29129487276077 -0.24794004857540, 
-4^727126628160E-02 . -0.38791200518608. 
0 194 8933987617 . -0.36395436525345. 
-0 40166890621 185 0.18578046560287. 
-0 46209335327148 -0.11633513122797. 
3"4i677090525627 0.29576727747917. 
-6.0622435063124EA)2 . 0.54043716192245. 
-0 43378528952599 . 0.37475159764290. 
-0 59756731987000 . -2909648977220 1E-02. 
-0*42078813910484 , -0.45302680134773. 
0'63260388374329 . 2.2896632552147E-02. 
0 47454643249512.-0.43297693133354. 
0 64523047208786 . 4.222498089075 1E-02. 
-041 148558259010 -0.49813264608383. 
S 6346767842770E-b2 . -0.63570129871368. 
0 52332174777985 , -0:35753899812698. 
-015361997485161 0.60483527183533. 
0:54942035675049 -0.27337777614594, 
0 24155624210835 , -0.55399334430695, 
O57553601264954 . -0.16218262910843. 
0 48504012823105 , 0.34690761566162. 

0.46577885746956 . -0.4002926647663 , 
-0.62349128723145 . -0.12451778352261. 
-0.30245357751846 . -0.59375917911530. 
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-0.64292871952057 , -0.2900027334 690 1, 
0.726071 1 1930847 , -0.19453555345535, 
0.65767973661423 , 0.46279323730060, 
0.85773074626923 .-7.9774774610996E-02, 
0.66653263568878 , 0.63700795173645, 
0.98371291160583 , 3.8463655859232E-02, 
-0.66836547851562 , -0.80670833587646, 
-1.0991182327271 , -0.15678122639656, 
0.96613508462906 , -0.66218781471252, 
-0.27184593677521 t 1.1993353366852, 

1.1098629236221 , -0.64718896150589, 

1.2801918983459 , 0.38049080967903, 

0.62276816368103 , 1.2329781055450, 
-0.47980448603630 , 1.3380939960480, 

0.58856463432312 , 1.3312312364578, 

0.56721973419189 , -1.3701412677765, 

-1.4011698961258 f 0.54447442293167, 

1.3742272853851 t 0.64058619737625, 

-1.4402492046356 , 0.49066120386124, 

1.3491046428680 , 0.69823592901230, 
-0.42755731940269 , -1.4469082355499, 
-0.73903095722198 » 1.2944310903549, 

-1.4206207990646 , 0.35585231252861, 

-1.2107812166214 , -0.76239341497421, 

-1.3619115352631 , 0.27647992968559, 

1.0996322631836 , 0.76832592487335, 
1.2723374366760 ,-0.19058661162853, 

0.96331912279129 , 0.75740247964859, 
1.1544399261475 , -9.9505998 19421 8E-02, 
-0.80496501922607 , -0.73074960708618, 

-1.0116646289825 , 4,7 18034993857 IE-03, 
-0.69000792503357 , 0.62838578224182, 
-0.84825158119202 , -9.2191 144824028E-02, 

0.63727533817291 t -0.43797433376312, 
-0.66910350322723 , -0.18957243859768, 

0.57503563165665 , -0.23856720328331, 

0.28575962781906 ( -0.47963011264801, 

3.5298123955727E-02 , 0.50607603788376, 

0.37911093235016 , -0.28557923436165, 

0.43339198827744 , 0.16655600070953, 

-9.2856124043465E-02 . -0.46805140376091, 
-0.36173328757286 , 0.36008602380753, 
•0.55U1116170883 t -9.266207367 1818E-02, 

034514288902283 . -0.28926151990891, 

0.26530939340591 ,-0.62695705890656, 

0.71201461553574 , -0.22391571104527, 
-0.69442093372345 , -0.41981008648872, 
-0.16683530807495 , -0.85803741216660, 
-0.55144131183624 , 0.75251996517181, 

0.12049946933985 , 0.97948336601257, 

0.80047225952148, 0.65618216991425, 

1.0733104944229 , -8.69883 14986229E-02, 

0.83770149946213 , 0.73083782196045, 
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-6.79078 1766 1762E-02 . -1. J 37 245 89347 84, 

-0.86384069919586 . -0.77313911914825. 
1.1698449850082 , -6.432849 16S7775E-02, 

-0.87872487306595 , -0.78181076049805. 
-1.1705212593079 , 7.6740458607674O-O2, 

-0.75661289691925 , 0.88237941265106, 

-0.10502639412880 . -1.1395549774170. 
0.87500441074371 , 0.69834393262863, 
1.0780709981918 , -0.14845751225948. 

-0.60881644487381 . 0.85695564746857. 
0.20570844411850 , 0.98799550533295, 

-0.49079838395119 . 0.82871961593628. 

-0.87198609113693 , 0.27489036321640. 
0.34792473912239 . -0.79089194536209. 
0.35360851883888 . 0.73334228992462. 

-0.18458327651024 , 0.74415004253387. 

-0.57590019702911 , 0.43903124332428. 
5.7776970788836E-03 . -0.68923103809357, 
OJ2798008918762 , 0.40390974283218, 

-0.62690776586533 . 0.18303082883358. 
0.22190481424332 , -0.61702972650528, 
0.55797028541565 . -0.37609174847603. 
0.70261883735657 , 3.4565746784210E-02, 

-0.56754517555237 , -0.48320794105530. 

-0.78116804361343 , 0.15341848134995. 

-0.40339776873589 , 0.75160259008408, 
0.84919923543930 . -0.33748552203178. 

-0.31929352879524 , 0.92272567749023, 

-0.51332873106003 , -0.90345317125320. 
0.23162299394608 , -1.0757941007614. 

-0.94100481271744 . 0.67701315879822, 
-1.2062630653381 , -0.14108607172966. 

-0.95936530828476 , 0.82507961988449, 
1.3102972507477 . 4.8358015716076E-02. 

-0.95462840795517 . -0.95657634735107, 
1.3848868608475 , -4.5906782 150269E-02, 
-1.0633902549744 , -0.93128252029419, 
1.4279329776764 , -0.14106513559818, 
1.1497660875320 , 0.88279163837433, 
0.23647691309452 , 1.4383053779602, 

-0.81110906600952 , U128539085388, 
0.33149531483650 , 1.4158697128296, 
-1.2524482011795 , -0.71695524454117, 
0.42545720934868 , 1.3614822626114, 

-0.60175400972366 , 1.2690199613571, 
1.2769542932510 , -0.51767551898956, 
0.46760299801826 , -1.2636756896973, 
0.60743457078934 , 1.1649851799011, 
1.2380973100662 , 0.31721931695938, 
0.69398546218872 , 1.0290712118149, 
0.15386630594730 , -1.1944670677185, 

-0.77655005455017 , -0.87338465452194, 
-1.8741229549050E-02 , -1.1353790760040, 
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-0.70263749361038 , 0.85432267 1S9026, 
-0.19654582440853 , -1.0637385845184, 
0.92648190259933 , 0.52192509174347. 
0.98265725374222 , -0.37526836991310, 
-0.99220246076584 , -0.33656212687492, 
0.55053806304932 , 0.89534473419189. 
1.0506724119186 , 0.15191125869751. 
0.71802806854248 , 0.80500024557114, 
1.1011143922806 , -2.6790 169999003E-02. 
0.87359136343002 , 0.71470791101456. 
-0.19459256529808 , -1.1428083181381. 
-1.0133891105652 , -0.62734061479568, 
-0.34699076414108 , -1.1751155853271, 
0.54547435045242 , -1.1340159177780, 
0.48006322979927 , 1.1975069046021, 
-1.2326081991196 , -0.47131305932999, 
0.59059274196625 , 1.2095856666565. 
0.40663030743599 , -1.3069387674332, 
0.67615872621536 , 1.21 1 1 126184464, 
1.3554904460907 , 0.35272744297981. 
0.73520946502686 . 1.2020268440247. 
-1.3775087594986 , -0.31040790677071, 
-1.1824613809586 , 0.76710045337677. 
-0.27996954321861 , 1.3730301856995. 
0.77210420370102 , 1.1527576446533. 
1.3428838253021 , 0.26121476292610, 
-0.75139236450195 . -1.1134694814682, 
1.2886705398560 . 0.25347581505775, 
-1.0653619766235 . 0.70698553323746. 
-1.2127152681351 . -0.25565749406815, 
1.0094038248062 , -0.64167791604996. 
1.1179939508438 , 0.26628923416138. 
0.94675332307816 , -0.55893445014954, 
0.28359553217888 , -1.0080425739288. 
-0.46277087926865 , -0.87873333692551, 
-0.88684558868408 , -0.30556473135948, 
-0.80680304765701 , 0.35761266946793, 
0.75870722532272 , 0.33003318309784, 
0.24814304709435 , 0.73252141475677, 
-0.62811541557312 , -0.35476726293564, 
0.65750586986542 , -0.13914376497269, 
-0.49959596991539 , -0.37754675745964, 
0.58338809013367 , -3.5333048552275E-02, 
0.37756574153900 , 0.39624589681625, 
0.51176452636719 , 5.8793801814318E-02, 
0.40890917181969 ,-0.26618802547455, 
-0.44414889812469 , -0.13911458849907, 
-0.41381889581680 , 0.16923630237579, 
-0.38192304968834 , -0.20201873779297. 
8.9968837797642E-02 , 0.40955415368080, 
-0.24452674388885 , 0.32629245519638, 
0.39503705501556 , -3.I020458787680E-02, 
0.26439002156258 , -0.27824369072914, 
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-0.36956810951233 , -5.6873057037592E-03, 
-0.26016449928284 , 0.23851068317890. 
-0.33284601569176 , - 1.90 10026007 89 1 E-02, 
-0.23125946521759 , 0.20754601061344, 
-0.28497639298439 , -8.6228372529 149E-03, 
0.17795602977276 , -0.18550240993500, 
-0.22646434605122 , 2.4967700242996E-02); 
double frl [5 14]; 

mainO 
( 

int crrorrc=0 f crrorim=0,numbcr l mf,inii l m l id[257] f i;[257] ( ij ( k; 
im jrcjim.l; 

double mctricp ) mctricm l minil,mini2jcy; 
double dcltphas, xrc, xim, curphas[129]; 

asmf.global declared"); 
asm ("declared:"); 

number =128; 

/************** rca d mc i & q components & find their amplitude, "frl" **+*+****+*/ 

for Q=0; j <= number, j++) ( 

curphasQ] = atan2(z[2*j+lU[2*jl); 

frl[2*j] = sqrt(z[2*j]*z[2*j].K42*j+l]*z[2*j+l));) 

asm(".global rcad.data"); 
asm ("rcad_data:"); 

/*********** Estimate die Group delay "frl" from die amplitude "frl" ************/ 

for (k= 1 ; k <= number ;k++) ( 

frl[2*k-l] = (frl[2*k] - frl[2*k-2])/frl[2*k]; 

frl[2*k-2]=0.0;) 
£rl[2*k] = 0.0; 

/************ estimate the group delay "frl" using the hilbcrt transform ************/ 
m=7; 

. asm('\global b_fft M ); 
asm ("bjft;"); 

fft(numbcr,m,frl); 
asm(".giobal ajft"); 
asm ("njfc"); 

for(k«l ; k<= numbcr/2; k++) [ 
jcy = frl[2*k-lj/numbcr, 
frl[2*k-l) a 0.0-frI (2*k]/numbcr, 
frl[2*kHcy;] 

for (b=numbcr/2+l; k <= number, k++) ( 
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icy = frl[2*k-l}/numbcr, 
frl{2*k-l]= frH2*k]/numbcr, 
frl[2*k]=0-jcy;] 

asm( n .ulobalbJfft"); 
asmC'bJfft: 1 '); 

ifftfoumbcr.nvfrl); 
asm('\giobai a_ifft"); 
asm( M aJfft:"); 



Thresholding & segmenting ***********************/ 

The purpose of this section is to determine the scgcmcnls of ^ received frequency frame 
wKTSdnulcd croup delay is larger than a certain threshold. These segments arc 
labelled 
****/ 

id[l]=l; 
m = 0; 
init = 0; 

for(k=l;k<=numbcr.k++)( 
if(( frl[2*k-l] > 0.07 )ll( rrl[2*k-l] < -0.07))( 
frl[k] = frl[2*k-l]; 
if(init~0)( 
m++; 
id[m] = k; 
inital;)) 
else ( 

frl[k] = 0.; 
if(init=l) ( 
m++; 

id[ni] = k; 

init = 0;)j 

if(inil— 1)( 
m++; 

id[m] =» number +1;J 

asmf.global threshold"); 
asmrihrcshold;"); 

y****»»«***************4i***Group delay Sign Algorithm********************* 4 "* 

The purpose of this part of ihc software is to determine the proper sign of estimated group 
delay. This is done through defining two Hamming distances one if we assucm the 
cstaimtcd group delay was positive and die other if it is assumed negative and choosing die 
sign that leads to die smaller distance 

*****************^^ 

for(k=l;k<=m;k++){ 
metricp ss 0.; 



12 



2064975 



APPENDIX C 

Copyright § 1991 AGT Limited 



mctricm = 0.; 

if( (id[kl=id[k+l]-l) I! (icltkl«*l[k+l]-2) )( 
s[k]=6; 
goto twenty;) 

for (j =id[k]+l; j<=id[k+l)-l; 

dcltphas — curphas[j] - curphas[j-l] - frl[j]; 

xrc » 1.4142136377*cos(ddtphas); 
xim = 1.4142136377*sin(dcltphas); 

mini! = l.-xrc; 

minil = minil*mimi; 
if(mini l>((l.+xrc)*(l .+xrc))) mini 1=(1 ,+xrc)*(l .+xrc); 

mini2 = l.-xim; 

mini2 = mini2*mini2; 
if(mini2>((l.+xim)*(l.+xim))) miiu2-(l.+xim)*(l.+xim); 

mclricp +<= minil + miiu2; 

dcliphas - curphas[j] - curphas[j- 1] + frlfj]; 

xrc » U4142136377*cos(dcllphas); 
xim « L4142136377*sin(dcllphas); 

minil = l.-xrc; 

minil « minil * minil; 
if(minil>((l.+xrc)*(L+xrc))) minil«(l.+xrc)*(l.+xrc); 

mini2 = l.-xim; 
mini2 = mini2 * mini2; 

if(mini2>((l.+xim)*(l.+xim))) mini2=(l.+xim)*(l.+xim); 

mctricm += minil + mini2; 

} 

s[k]=M; 

if(mctricp < mctricm) s[k]=l; 
. /* printf( M s[%d]=.%d\n l, ( k,s[k]);*/ 
twenty: k++;) 

for (k»l ; k<=m; k++) { 
for 0 *id[kj; j<=nd[k+l]-l; j++) frllj] *» s(kj; 
/*k++;*/ ) 

asmr.globalajs"); 
asm^aJs:"); 

Decode the Signal ***+*************♦**********/ 
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1 = 1; 

while (1 < numbcr+l)( 

dchphas = curphasll] - curphas[i-l] - frl[l]; 

xrc = cos(dcltphas); 
xim = sin(dcliphas); 

jrc = xrt/fabs(xrc); 
jim = xim/fabs(xim); 
if (jrc=jim) (jrc = 0 -jrc; 
jim = 0 -jim;) 

if(l>Dl 
if(jrc !=jj2[2*l]) crrorrc++; 

if(jim != jj2[2*l+l]) crrorinvH-; 
] 1+-K 

printf("%d %d %cNi , \i l crrorrc ( crrorini); 

asmCglobalcnd"); 
asm("cnd:"); 

} /*cndof main*/ 
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40 

THE EMBODIMENTS OF THE INVENTION IN WHICH AN EXCLUSIVE 
PROPERTY OR PRIVILEGE IS CLAIMED ARE DEFINED AS 
FOLLOWS : 

1. a transceiver including a transmitter for 

transmitting electromagnetic signals and a receiver 
for receiving electromagnetic signals having amplitide 
and phase differential characteristics, the 
transmitter comprising: 

an encoder for encoding information; 

a wideband frequency division multiplexer 
for multiplexing the encoded information onto wideband 
frequency channels; 

a processor having an output for applying a 
Fourier transform to the multiplexed information to 
bring the information into the time domain for 
transmission; 

a low pass filter on the output from the 

processor; 

a local oscillator for upconverting the 
multiplexed information for transmission; and 

further including, in the receiver of the 
transceiver; 

a bandpass filter for filtering the received 
electromagnetic signals; 

a local oscillator for downconverting the 
received electromagnetic signals to produce output; 

a sampler for sampling the output of the 
local oscillator to produce sampled signals to the 
channel estimator; 

a channel estimator for estimating one or 
both of the amplitude and the phase differential of 
the received signals to produce as output one or both 
of an estimated amplitude and an estimated phase 
differential respectively; and 
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a decoder for producing signals from the 
sampled signals and the output from the channel 
estimator. 

2. The transceiver of claim 1 further 
including, in the receiver of the transceiver: 

a deprocessor for applying an inverse 
Fourier transform to the samples output from the 
sampler. 

3. The transceiver of claim 1, further 
including, in the receiver of the transceiver: 

a power controller before the sampler for 
monitoring the power of the received signal and for 
controlling the power of the signal. 

4. The transceiver of claim 1, further 
including, in the transmitter of the transceiver: 

a pre-distorter before the processor, the 
pre-distorter being connected to the channel 
estimator, for predi storting a signal to be 
transmitted with one or both of the estimated 
amplitude or the estimated phase differential. 

5. The transceiver of claim 4 in which the 
power controller is also connected to the pre- 
distorter for controlling the power of the signal to 
be transmitted. 

6. The transceiver of claim 1 further 
including: 

means to modify the received signal with one 
or both of the estimated amplitude and phase 
differential respectively. 
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7. A method for allowing a number of wireless 

transceivers to exchange frames of information, the 
method comprising the steps of: 

multiplexing a first frame of information 
5 over a number of frequencies within a frequency band 

at a first transceiver to produce multiplexed 
information; 

processing the multiplexed information at 
the first transceiver; 
10 transmitting the processed information to a 

second transceiver using a carrier frequency f c ; 

receiving the processed information at the 
second transceiver; 

processing the processed information at the 
15 second transceiver during a first time interval; and 

the width of the frequency band being chosen 
so that neither carrier nor clock recovery is required 
at the second transceiver. 

20 8. The method of claim 7 in which the frequency 

band is formed from a first set of K x points and a 
pair of tail slots each having K 2 points, each of the 
points being separated by a frequency range Af, the 
second transceiver has a maximum expected clock error 

25 x T / where T is the duration of one time domain sample, 

the information is multiplexed over a number M of 
levels, and K x is selected such that 2it%/K x < ti/M. 

9. The method of claim 8 further including 

30 transmitting a second frame of information from the 

second transceiver to the first transceiver within the 
same frequency band. 
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10. The method of claim 8 in which K 2 is 
selected so that the out of band signal is less than 
a given level. 

11. The method of claim 8 in which the first and 
second transceivers have an expected maximum relative 
velocity V, the first and second transceivers have 
carrier frequencies with a frequency offset from each 
other of f Q , the carrier frequency has a corresponding 
travelling wavelength k and Af is selected so that 
[V/(XAf) + f Q /Af] is less than or equal to a 
preselected sampling error. 

12 . The method of claim 8 in which processing 
the mutliplexed information at the second transceiver 
further includes calculating the mean of the phase 
shift due to sampling error by summing the estimated 
phase differential of the received signal. 

13 . The method of claim 12 in which the mean of 
the phase shift due to sampling error is divided by K x 
and the result removed from the phase differential of 
the received signal. 

14. A transmitter for transmitting 
electromagnetic signals having amplitide and phase 
differential characteristics, the transmitter 
comprising: 

an encoder for encoding information; 

a wideband frequency division multiplexer 
for multiplexing the encoded information onto wideband 
frequency channels; 

a processor having an output for applying a 
Fourier transform to the multiplexed information to 
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bring the information into the time domain for 
transmission; 

means for operating on the multiplexed 
information to introduce an arbitrary phase into the 
multiplexed information; 

a low pass filter on the output from the 

processor ; and 

a local oscillator for upconverting the 
multiplexed information for transmission. 

15. The transmitter of claim 14 further in 
combination with a receiver for receiving 
electromagnetic signals, the receiver comprising: 

a bandpass filter for filtering the received 
electromagnetic signals; 

a local oscillator for downconverting the 
received electromagnetic signals to produce output ; 

a sampler for sampling the output of the 
local oscillator to produce sampled signals to the 
channel estimator; 

a channel estimator for estimating one or 
both of the amplitude and the phase differential of 
the received signals to produce as output one or both 
of an estimated amplitude and an estimated phase 
differential respectively; and 

a decoder for producing signals from the 
sampled signals and the output from the channel 
estimator. 

16. The transmitter and receiver combination of 
claim 15 further including, in the receiver: 

a deprocessor for applying an inverse 
Fourier transform to the samples output from the 
sampler. 
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17. The transmitter and receiver combination of 
claim 16, further including, in the receiver: 

a power controller before the sampler for 
monitoring the power of the received signal and for 
controlling the power of the signal. 

18. The transmitter and receiver combination of 
claim 15, further comprising: 

a pre-distorter before the processor, the 
pre-distorter being connected to the channel 
estimator, for predistorting a signal to be 
transmitted with one or both of the estimated 
amplitude or the estimated phase differential. 

19. The transcmitter and receiver combination of 
claim 18 in which the power controller is also 
connected to the pre-distorter for controlling the 
power of the signal to be transmitted. 

20. The transmitter and receiver combination of 
claim 15 further including: 

means to modify the received signal with one 
or both of the estimated amplitude and phase 
differential respectively. 

21. A method for allowing a number of wireless 
transceivers to exchange frames of information 
symbols, the method comprising the steps of: 

multiplexing a first frame of information 
symbols over a number of frequencies within a 
frequency band at a first transceiver to produce 
multiplexed information; 

processing the multiplexed information at 
the first transceiver, including introducing an 
arbitrary phase into the information symbols; 
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transmitting the processed information to a 
second Transceiver using a carrier frequency f c ; 

receiving the processed information at the 
second transceiver; and 
5 processing the processed information at the 

second transceiver during a first time interval , 

22. The method ot claim 21 in which the width of 
the frequency band being chosen so that neither 

10 carrier nor clock recovery is required at the second 

transceiver* 

23. The method of claim 22 in which the 
frequency band is formed from a first set of K x points 

15 and a pair of tail slots each having points, each 

of the points being separated by a frequency range Af , 
the second transceiver has a maximum expected clock 
error %T, where T is the duration of one time domain 
sample, the information is multiplexed over a number 

20 M of levels, and K x is selected such that Zrtx/K^ < 

tt/M. 

24. The method of claim 23 further including 
transmitting a second frame of information from the 

25 second transceiver to the first transceiver within the 

same frequency band. 

25. The method of claim 23 in which is 
selected so that the out of band signal is less than 

30 a given level. 

26. The method of claim 23 in which the first 
and second transceivers have an expected maximum 
relative velocity V, the first and second transceivers 

35 have carrier frequencies with a frequency offset from 



JUL-06-1998 13:49 FROM THOMPSON LAMBERT 



TO 



18197709265 P. 03 



47 20 6 4 9 7 5* 

each other of f Q , the carrier frequency has a 
corresponding travelling wavelength X and Af is 
selected so that [v/UAf) + f 0 /Af] is less than or 
equal to a preselected sampling error. 

S 

27. The method of claim 23 in which processing 
the mutliplexed information at the second transceiver 
further includes calculating the mean of the phase 
shift due to sampling error by summing the estimated 

10 phase differential of the received signal. 

28. The method of claim 27 in which the mean of 
the phase shift due to sampling error is divided by K x 
and the result removed trom the phase differential of 

15 the received signal . 

29. A transmitter for transmitting 
electromagnetic signals having amplitide and phase 
differential characteristics, the transmitter 

20 comprising; 

an encoder for encoding information; 
a wideband frequency division multiplexer 
for multiplexing the encoded information onto wideband 
frequency channels; 
25 a processor having an output for applying a 

Fourier transform to the multiplexed information to 
bring the information into the time domain for 
transmission; 

a low pass filter on the output from the 
30 processor; and 

a local oscillator for upconverting the 
multiplexed information for transmission. 
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